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SECTION 1 
(This Section Is Unclassified) 


INTRODUCTION 
1.1 GENERAL 


In a previous effort, we evaluated dynamic pressure impulse, Iq, for 
ideal /near-ideal blast waves by showing a correlation with wheeled vehicle 
(1/4 ton trucks) displacement (Reference 1). We found that there is a one- 
to-one correspondence between Iq and vehicle displacement and concluded that 
from known vehicle displacements the dynamic pressure impulse can be esti- 
mated with reasonable accuracy. 

Here, we present the results for evaluating Iq, not only for ideal/near- 
ideal blast waves but also for Tq for non-ideal blast waves. We combined the 
data from our first evaluation with 1/4 ton truck exposures to non-ideal 
blast waves along with that of 2-1/2 ton truck exposures to both non-ideal and 
ideal/near-ideal blast environments. 

Additionally, we evaluated damage to wheeled vehicles as a function of 
Vehicle displacement. For a given displacement, the damage can then be corre- 


by simply using Sachs scaling. The damage itself as a function of range, how- 
ever, cannot be scaled by use of a constant exponent of the weapon yield such 
as yo.4, The exponent of W will vary depending on the burst height of the 
weapon. 

To keep this report from being too cumbersome we do not include the 
damage versus displacement here. Instead another report was prepared to 
Present this information. We also feel that in this way the report on Iq will 
be of special interest to groups dealing with airblast phenomena while the 
report on damage will be of special interest to groups dealing with targeting 
and vulnerability/survivability. This report on Ig presents the procedures 
for correlating Iq versus displacement, the effect of surface conditions on 
Iq, the construction of Ig versus range curves as well as the construction of 
HOB charts for iso-I,. 


1.2 OBJECTIVES 

Our main objective was the evaluation of dynamic pressure impulse for 
ideal/near-ideal and non-ideal blast waves using wheeled vehicle response 
along with gage measurements. Other objectives were to correlate damage 
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with vehicle response, establish the influence of diffraction flow on dis- 
placement and evaluate the peak overpressure effects on damage to wheeled 
vehicles. 


1.3 BACKGROUND 

The exposure of wheeled vehicles (drag-type targets) on past nuclear and 
HE tests served several purposes such as troop indoctrination, damage evalua- 
tion and verification of damage predictions. Different types of vehicles were 
exposed under a variety of conditions. The vehicles included 1/4 ton, 3/4 ton, 
2-1/2 ton, and 5 ton trucks. The majority of vehicles exposed were 1/4 ton 
trucks and on several occasions they were used as gages to evaluate the rela- 
tive effects of nuclear blast on drag-type targets. The information on 
wheeled vehicle exposures has been compiled and reevaluated. A report has 
been prepared which covers the exposure conditions, vehicle damage, and dis- 
placements (Reference 2). 

Damage to and displacements of drag-type targets under some nuclear burst 
conditions was far in excess of damage and displacements for other nuclear 
burst conditions. This was dramatically demonstrated on events of Encore-9 
and Grable-10, Upshot/Knothole. The maximum damage of vehicles (1/4 ton 
trucks) exposed on Encore-9.was M.6 (Moderate II) and the maximum displace- 
ment was about 11 metres. On Grable-10, however, the 1/4 ton trucks were 
completely dismembered and displacement of parts (engine, chassis, wheels) 
was greater than 300 metres. The exposure of vehicles on these two events 
was based on peux overpressure values and it was assumed that the other blast 
parameters could be calculated by using the Rankine-Hugoniot relations. But 
the pressure waveforms measured on Encore-9 (Mach Region) were primarily 
steep-rising shock fronts (near-ideal), while the pressure waveforms on 
Grable-10 were somewhat distorted showing rounded fronts with high frequency 
oscillatiuns and in some cases secondary seaks were observed (non-ideal). 
There was ns corselation between overpressure and dynamic pressure on Grable- 
10. Thic difference in damage and displa:ement was then attributed to the 
precursur dust-laden blast wave formed on Grabie-10. The dust momentum 
increased the total loads on the vehicles. <“imilar results were obtained 
on ether precursor-forming events. Tius one of the criteria used for dis- 


tinauishlng between ideal,/near-icseal and non-ideal blast environments was the 
waveform — steep rising shock front for ideal/near-ideal and a rounded front 
with secondary reaks for non-ideal. The non-ideai could be further divided 
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this was that dynamic pressures were greater on heavy dust events than on 
light dust environments. 

When vehicles were used as gages they were located adjacent to blast 
line instrumentation which consisted of overpressure, pitot-tube, differential 
pressure, total head, GREG and SNOB gages. In a precursor dust-laden blast 
environment the dust registry coefficient was known for the GREG and SNOB 
gages only. Many of the measurements for deriving dynamic pressure versus 
time in this sort of an environment were suspect. A procedure was developed 
to try to retrieve the suspect data (Reference 3). The technique for retriev- 
ing the suspect data was based on the comparison of measured data from two 
gages — one gage with a known registry coefficient and the other with an 
unknown coefficient. This procedure was not successful since there were a 
limited number of measurements with gages having a known registry coefficient. 
Wheeled vehicles on the other hand responded by translational and rotational 
motion to the actual blast load environments. The displacement of the vehicle 
should, therefore, be a measure of the blast loads acting on the target. 

Past studies have resulted in calculational capability to compute 
rotational and translational motion of drag-type targets (References 4 and 5). 
The technique of Reference 4 calculates either translational motion or rota- 
tional motion. The technique of Reference 5 which is a multidegree-of-freedom 
code calculates both translation and rotation. However, for both techniques 
the calculation for rotation is limite up to the point of overturning. The 
important blast parameter considered ° . both techniques is dynamic pressure 
impulse. For fractional KT weapons, diffraction loading due to overpressure 
influences the motion. These calculations provide an insight for proper 
interpretation of the vehicle response to the blast loads. 

In our initial evaluation, we were concerned about establishing the 
values of dynamic pressure impulse for the ideal/near-ideal blast waves to 
be correlated with vehicle displacements. Two approaches were used to 
establish these values. One approach was based on averaging the measured 
overpressure and positive phase duration as a function of range and then 
combining this with the Rankine-Hugoniot relations together with empirically 
derived equations to compute the peak dynamic pressures and the peak dynamic 
pressure impulse, respectively. Another approach was the use of theoretical 
and analytical computations for static pressure impulse, Ip, and Iq combined 
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with the average values of measured peak overpressure, P,, and static pressure 
impulse. The tatter approach, which was the method selected, permitted the 
construction of a set of curves showing the ratio of Ip/Iq versus P, for 
several burst heights. The I, was determined using the average measured values 
of ip and Ps on each event. Details of this approach are given in Appendix A. 
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SECTION 2 
(U) BLAST MEASUREMENTS AND VCH™C'.E EXPOSURES - NUCLEAR AND HE EVENTS 


2.1 (U) TECHNICAL DATA BASE 


2.1..  (U) Opera. ions/Events 

(U) Over 400 vehicles wert: exposeu on 19 nuclear and HE events. Blast 
measurements of one type or anothe: were also obtained from these events. Of 
the more than 400 vehicles exposed, only 167 were used for correlation with 
I, The vehicles excluded for this correlation were those located on the 
asphalt and water surfaces and those located in the Desert Rock sector, and 
also, vehicles which were constrained to rotate without sliding. Blast over- 
pressure measurements were obtained on the majority of these events. Only on 
two events, Fncore-9 and Wasp-1, were the dynamic pressure measurements un- 
ava7lable. 

(U) The listing of the Operation/Events are given in Table 2.1. The 
information provided includes the yield, actual burst height, scaled burst 
height, scaling factors, and ambient environment. Ten of these events are 
‘additions to the events used in our first evaluation. These added events all : 
took place at Nevada Test Site; most produced strong precursor actions (non- 
ideal blast waves). The nomenclature given for each event refers to the code 
name and to the shot number for that Operation. 


2.1.2 (U) Blast Wave and Vehicle Displacement Collations 
(U) The nuclear and HE events, designated as ideal/near-ideal blast 
waves were the assigned nine events (first evaiuation) and the additional two 
events: Wasp-] and Wasp Prime-9. Our initial guidance for cutegorizing the 
precursor nuclear events was based on the information given in Reference 6, 
j.e., whether they were light dust-laden or heavy dust-laden blast waves. 
The sorting of events into one of the three categories is as follows: 


MN Ideal/Near-Ideal Blast Waves 


Upshot/Knothole: Encore-9 


i! 
Castle: Kocen-3 and Nectar-6 
Teapot: Wasp-1 and Wasp Prime-9 
Redwing: Lacrosse-1, Zuni-3 and Yuma-4 


HE: Canadian 20-Ton, 100-Ton and Dice Throw 
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Light Dust-Laden Blast Waves 
Teapot: Moth-2, Turk-4, Bee-6, and Apple I-8 


Heavy Dust-Laden Blast Waves 

Upshot/Knothole: Grable-10 

Teapot: Met-12 and Apple II-13 

Plumbbob, Smoky-15 
Some of the events within the light dust and heavy dust categories, however, 
were interchanged as will be described later in this report. The blast mea- 
surements of total pressure and overpressure as a function of time were also 
gatnered and categorized accordingly. 
(U) The wheeled vehicle exposure data were sorted into these three 
categories. The vehicle data pertinent to this evaluation were the location, 
orientations, displacements and damage. Vehicles were further sorted according 
to surface conditions, i.e., desert, water, asphalt. 


2.2 (U) ANALYTICAL APPROACH 


2.2.1 (U) Blast Data Reduction Procedures 
(U) We have described the procedures used for deriving Ig for ideal/ 
near-ideal blast conditions (see Background and Appendix A). In order to » 
obtain the initial best estimates of Ig for the non-ideal blast conditions we 
reduced the data independently of the previous data reduction procedures. The 
procedure we followed was to integrate the total pressure-time and the static 
pressure-time waveforms to obtain the total pressure impulse and the static 
pressure impulse. A check was inade for certain intervals of the total pressure- 
time and the static pressure-time curves to determine the range of Mach numbers 
for those intervals. Where appropriate the Mach compressibility corrections 
were made to the difference between the total pressure impulse and static 
pressure impulse for those intervals. The sum over these intervals yielded 
the dynamic pressure impulse. The details of this procedure are presented in 


Se 


-~ 
— 


oe oe 
— 


fs 


- 
. 


Appendix B. 

(U) We made checks of this procedure by comparing the Iq values we 
derived with those reported by BRL on Plumbbub, Smoky-15 (Reference 7) and on 
Dice Throw (Reference 8). The Macn numbers ranged as high as 2.2 in these 
comparisons. Waveforms on Smoky~15 displayed a considerable amount of oscil- 
lations while the waveforms selected from Dic~ Throw were reasonably smooth. 
On Smoky-15, at the first stations, we assu’ that the gage was plugged since 
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curve corresponding to the end of the positive-phase on the static pressure- 
time curve. We assumed that the error contribution to the total pressure 
impulse would be small at the tail-end of the waveform. The comparisons of 
Iq are shown in Table 2.2. The maximum difference between these two values 
is about 15 percent. 


2.2.2 (U) Least Squares Fit Procedures 
(U) All of the problems evaluated where least squares fits were applied 
made use of the TI-59 Programmable Calculator. The problems evaluated were: 


UNCLASSIFIED 
the total pressure-time curve did not return to the base line. For computing 
the total pressure impulse we terminated the integration of the pressure-time | 
1) Dynamic Pressure Impulse vs. Displacement [subcases included 
WWII and M38A1 (1/4 ton trucks) and the M35 and M135 (2-1/2 ton trucks) in 
both side-on and face-on orientations]. 
2) Damage vs. Displacement (several subcases). 
3) Displacement vs. Scaled Ground Range (subcases included 
side-on and face-on for 1/4 ton trucks with various height of burst conditions 
including light and heavy dust). ; 
° 4) Dynamic Pressure Impulse vs. Scaled Ground Range (subcases 
included air burst and surface burst as well as light and heavy dust condi- 
tions). 
5) Scaled Displacement vs. Scaled Ground Range (WWII, 1/4 ton 
trucks, side-on and face-on, surface burst conditions). 
6) Scaled Dynamic Pressure Impulse vs. Scaled Ground Range 
(subcases included air burst and surface burst as well as light and heavy 
dust conditions). | 
The functional forms used for least squares fits for these 
preliminary evaluations were as follows: | 


PROBLEM 1: Iq = BO” or D = Aly” 
where: a and b are coefficients 
m and n are exponents 
D is displacement (metres) 


a ot oa 2 
PROBLEM 2: da ag * a,D + aD » D < 0, 


da = 10, D>D, 
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damage 

displacement 

value of displacement for which curve fit 
(first) reaches line da = 10 


where: da 
D 


PROBLEMS 3 through 6: £n y = ay + a,x + ayx” 


where: y = ordinate, Displacement or I 
x = independent variable, Scaled Ground 
Range 
(U) A detailed description of the least squares fit procedures is 


given in Appendix C. 

As stated previously, the above functional forms were used for 
the preliminary evaluations. Later, we used other functional forms and these 
will be defined in other Sections of this report. 


2.2.3 (U) Displacement - Dynamic Pressure Impulse Least Squares Fits 
(U) In Appendix A we have described how we obtained values of Iq for 
the ideal shots used in this report. In Appendix B we have described in con- 
siderable detail how we obtained values of Ig from gage data for the non-ideal 
shots of present concern. The results are contained in Table 2.3. Figures 
2.1 through 2.4 show plots of these values of Iq versus vehicle displacement 
for the WWII and M38A1 one-quarter ‘ton trucks in side-on and face-on orienta- 
tions. * 
(U) We obtained least squares fits for both Ig versus displacement and 
for disnlacement versus Ig in each case following the procedures described in 
Appendix C. The results are given in Table 2.4. Similar results are given for 
the REO M35 2-1/2 ton trucks in Figures 2.5 and 2.6 and Table 2.4. (For the 
REO M35 improved or "first iteration” values of I, were used as described in 
Section 2.2.4.) 
(U) We will discuss these results in detail in Section 3.1. Here we 
note, however, that while the data exhibit considerable scatter, the curve 
fits to the data are well determined. Thus, on the average, Iq can be esti- 
mated with fair accuracy from a measurement of vehicle displacement — the 
fitted curve serving as a calibration and the vehicle serving as a gage. The 
results* in Table 2.4 show that values of Iq inferred from displacement mea- 
surements are of comparable accuracy with the gage-reduced data for the non- . 
ideal shots discussed in Appendix B. 


(U) *Note the listed: values of ET: 
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2.2.4 (U) Improved or "Iterated" Values of I, 
(U) One of the objectives of this report is to obtain the best, most 
reliable possible values of dynamic pressure impulse for non-ideal nuclear 
shots. Since vehicle displacement is independent of the detailed cause of 
the impulse to which the vehicle is subjected, i.e., whether it is caused 
solely by airblast or by an air-dust mixture, and since displacement-inferred 
values of impulse appear to be comparable in accuracy with gage-inferred 
values, it appears logical to use displacement-inferred results to improve 
our estimates of dynamic pressure impulse. Therefore, for each instance in 
this report where we have one or more 1/4 ton trucks exposed on a non-ideal 
shot, side-on and/or face-on, at the same ground range as that of a gage loca- 
tion — in most cases gages and vehicles were placed at the same ground range 
— we have averaged the values of dynamic pressure impulse inferred from the 
measured displacement by use of the dynamic pressure impulse versus displace- 
ment data fits. (Note that these fits are based on both ideal/near-ideal and 
non-ideal nuclear shot data, and include some HE shot data also.) This pro- 
vides, then, a single, mean value of dynamic pressure impulse inferred from 
vehicle displacements. 
(U) Finally, since the gage-inferred and vehicle displacement-inferred 
values of dynamic pressure impulse appear to be of comparable accuracy, we 
have averaged the (average) gage-inferred value with the (average) displace- 
ment-inferred value to obtain a "final" average or "grand" average value of 
dynamic pressure impulse. We note that in most cases the number of gage- 
inferred values producing the gage average result is 2 or 3 (i.e., 2 or 3 
of the 4 possible values acceptable); the displacement-inferred result is 
usually based on the average of 2 results, one side-on and one face-on. 
However, the displacement results are independent whereas the gage-inferred 
results are partially dependent in a complicated way, as we have noted. 
(U) The improved estimates of dynamic pressure impulse obtained for 
non-ideal shots by this procedure are given in Table 2.3. These results, 
along with our original estimates for Ig for ideal and near-ideal shots, are 
used in Sections 3 and 4 of this report. 
(U) Having obtained the above results for .1e non-ideal shots as 
described, we could consider replotting displacement versus Ig using the 
improved Iq results with the expectation that data scatter would be lessened 
(since the average curve for the originally plotted results is used to improve 
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the estimates of some of the points now being replotted)! However, the dimi- 
nution in the data scatter would actually be quite small: the scatter of the 
ideal/near-ideal data points about the curves is as great as for th- non-ideal 
(and there are more ideal/near-ideal data points). This may seem *“‘ogical: 
gage data should be better for ideal shots since the gages are not si bject to 
clogging by dust and pressure-time waveforms are freer of irregularities. It 
follows from this that the ideal gage data should also be better than vehicle 
displacement-inferred results since the latter were found to be comparable in 
accuracy to the gage-inferred values of Iq for non-ideal shots. For the fact 
that the ideal nuclear shot results show as much data scatter in the displace- 
ment-Iq plots as do the non-ideal we offer the following explanation: 

1) there is a greater variety of ideal/near-ideal shots, 
especially with regard to the range of weapon yields in the data base; 

2) our calculation of four gage-derived values of Ig, largely 
independent from the standpoint of data reduction procedure (see Appendix B), 
and intercomparison of the calculated values, after culling obviously bad 
gage data, resulted in elimination of most of the poor data (whether caused 
by a gage result or the data reduction process) before averaging to obtain 
the accepted gage-derived result; (the latter was then averaged with the 
average displacement-derived result to obtain our improved estimate). 


Whether or not the above explanation accounts entirely for it, the fact is 

that ideal and non-ideal data have about the same scatter on the displacement 
versus Iq plots. Thus, while initially it seemed reasonable to improve our 
estimates of Ig for the non-ideal shots by using the vehicle displacement data 
— the "calibration" of the vehicle as a gage being determined by both the 
ideal/near-ideal data along with our initial estimates of non-ideal data (gage 
-inferred) — it now appears equally reasonable to use the vehicle displacement 
data to improve Ig estimates for the ideal/near-ideal shots: the results for 
the latter are apparently not as accurate as might have been supposed. The 
same rationale justifies the procedure: the central curve in the displacement 
“Iq plots is well determined; the data point scatter about this curve is fairly 
large but quite uniform. So while displacement-inferred values of Iq have 
considerable scatter (there being several uncontrolled variables which affect 
displacement, (see Section 3.1), their accuracy is comparable with that of the 
gage-inferred and/or calculated results for the ideal/near-ideal shots. 
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values of I_ to improve our estimates in order to obtain the best possible 
estimates of I, for construction of scaled HOB contours. To distinguish 
these results from the previous results, we now refer to the first such 
results (Table 2.3), applicable only to non~ideal shots, as "first iteration" 
results. Our second use of this procedure, which applies to both idea? and 
non-ideal shots and which starts with the best previous estimate (i.e., initial 
estimate for ideal/near-ideal shots, “improved” or "first iteration" results 
for non-ideal shots) — we now refer to as the "second iteration". We note 
here the changes produced in this iteration were quite modest. However, as 
will be seen in Section 4.1, it is quite difficult to determine acct” itely 
the shapes of scaled HOB I, contours so that the best possible data should 

be used for this purpose. Section 4.1 is the only part of this report in 
which we use "second iteration" values °f I.. The second iteration i.e., 


final best estimates of I. are given in Tables 4.1 through 4.4. 


(U) -oliowing the above reasoning we again used displacement-inferred 
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TABLE 2.2 (U) Comparison of Dynamic Pressure Impulse 
UNCLASSIFIED 


Plumbbob, Smoky-15 


baleutasee BRL (Report WT-1407) 
I 
RANGE 4 % 
Metres psi-sec kPa-sec psi-sec kPa-sec 
457 20.02 138 21.99 152 
789 3.44 23.7 3.85 26.5 
897 2.29 15.8 2.02 13.9 
1038 2.28 15.7 2.26 15.6 
Dice Throw 2 
Calculated BRL (Report 0NA4377P-1) 
I I 
RANGE %0 90 
(Metres ) psi-sec kPa-sec psi-sec kPa-sec 
134 1.56 10.79 1.45 9.98 
195 0.59 4.07 0.53 3.67 
375 0.096 0.66 0.082 0.57 
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TABLE 2.3 (U) Dynaric Pressure Impulse vs. Range, 
Non-Ideal Blast Waves 


UNCLASSIFIED 
Displacement "Grand" Avg 
Ground Gage Reduced Inferred (First Iteration 
Range Result (avg) I, (avg) Result) 
Event (mn) (kPa-sec) (kPh-sec) (kPa-sec) 
Grable-10 736 5.1 5.16 5.13 
844 2.5 2.87 2.69 
Moth-2 287 6.75 -- 6.75 
366 3.35 -- 3.35 
41] 2.54 2.31 2.42 
457 1.88 1.83 1.85 
518 1.35 -- 1.35 
Turk-4 716 20.7 -- 20.7 
914 11.8 -- 11.8 
1030 7.95 7.42 7.69 
1128 5.50 5.82 5.66 
1280 2.32 -- 2.32 
Bee-6 ‘ 549 6.1 10.1 8.10 
. 610 3.75 4.59 4.17 
701 2.25 2.40 » 2.32 
777 1.35 2.47 1.91 
Apple I-8 527 17.3 -- 17.3 
622 9.10 -- 9.10 
902 1.46 2.44 1.95 
991 0.96 1.44 1.20 
1128 0.41 0.338 0.38 
Met-12 610 48.0 38.0 43.0 
686 37.2 38.0 37.6 
762 17.2 16.3 16.8 
838 41.5 16.1 13.8 
914 3.4 3.62 3.51 
Apple II-13 518 41.9 -- 41.9 
625 34.0 -< 34.0 
808 19.3 -- 19.3 
914 12.1 8.54 10.3 
1006 7.4 3.68 5.54 | 
1128 4.8 3.15 3.98 
1219 2.48 2.62 2.55 
smoky-15 726 22.3 34.7 28.5 
84) 14.7 17.3 16.0 
897 18.2 6.81 12.5 
1038 5.96 7.84 6.9 
1181 3.96 3.54 3.75 
14 
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(U) Dynamic Pressure Impulse versus Displacement - 1/4 Ton 


Trucks, WWII, Side-on, Desert Data Curves, Non-Desert 


Data Points 


FIGURE 2.1-1 
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FIGURE 2.2 (U) Dynamic Pressure Impulse versus Displacement - 1/4 


Ton Truck, 438A1, Side-on 
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3.1 (U) ANALYSIS OF VEHICLE MOTICN 


3.1.1 (U) Dynamic Pressure Impulse versus Displacement Results 
(U) The results for Dynamic Pressure Impulse versus Displacement are 
given in Figures 2.1 through 2.6 and Table 2.4. 
(U) The log-log plots indicate a rather uniform scatter of the data 
about the straight line fits. The values of r, the coefficient of correlation, 
are given in the tables. While r is a measure of the goodness of fit of the 
selected fitting function with respect to the data (see Appendix C, "Least 
Squares Fitting Procedures"), it appears that better values would not be expec- 
ted with a different fitting function; the data scatter preclude that. 
(U) There is more data for the WWII than for the M38Al vehicle. The 
fact that there is more scatter in the WWII data stems from the fact that a 
largernumber and variety of shots give rise to the data for tne WWII. Also 
for some ov the greater displacements, large resultant errors are introduced 
because of the dismemberment of the vehicles, particularly for the side-on 
orientations. There is some tendency for points from given shots to be above 
or below the fitted curve. Since this tendency has superimposed upon it the 
data point scatter associated with each shot, the total scatter is greater 
than for the M38Al vehicle. Thus, for the WWII, side-on plot, Figure 2.1, the 
data corresponding to shots Koon-3, Nectar-6, Canadian 20-Ton and 100-Ton have 
somewhat smaller displacements for a given value of Ig than the curve fit while 
data from shots Wasp Prime-9, Yuma-4, Bee-6 and Apple I-8 have somewhat higher 
displacements. Similarly for the WHII, face-on plot, Figure 2.3, (except that 
here no Canadian 20-Ton and 100-Ton points are shown — see next paragraph). 
Thus, a small HOB effect is present in the data. The data points for the M38Al, 
with very limited shot variety, result in a rather uniform scatter of the data 
for individual shots about the curve fit as best we can judge this from the 
small number of data points. i 
3.1.1.1 (U) Effect of Diffraction 
In the side-on orientation both vehicles usually turned 
over, thus causing displacements, as measured — before and after position of 
vehicle canter of mass — to be larger for side-on than for face-on orienta- 
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tion, except possibly for very large displacements where they appear to be 
equal for the WWII vehicle. In the face-on orientation the vehicle turned over 
infrequently in the case of the WWII and not at all for the M38Al. A few 
points corresponding to low vield shots are not shown on the plots nor are 
they included in the curve fits. These show the effects of diffraction. Gata 
for very low yield shots (10 tons) exhibit much higher values of visplacement 
for small values of dynamic pressure impulse than the plotted data. This is 
because the blast wave is so short that the vehicle received almost al] of the 
dynamic pressure impulse on its struck side while also subjected to a "Static" 
overpressure before any significant counterpressure has affected the opposite 
side. In the face-on data this effect is evident also for the Canadian 20-Ton, 
100-Ton and Yuma-4 (fractional KT, see Table 2.1) shots. The effect for these 
shots is not evident in the side-on data. We believe the reason for this is 
that in the side-on case vehicle turnover caused comparable displacements for 
Tow values of Ig for the higher yield shots, thus concealing the diffraction 
effect. 
3.1.1.2 (U) Partial Interpretation of Results 
(U) When a vehicle is struck by a blast wave, it will acquire 
a velocity proportional to the impulse received (if we ignore -the effect of 
the motion-resisting force during the time of delivery of the impulse). -There- 
fore, we expected to find a good correlation between displdcement and Iq. In 
later discussion we will note that damage correlates with displacement wehreas 
in a wide variety of studies involving blast damage to equipment, damage seems 
to depend upon both pressure and impulse. We, therefore, made plots of dis- 
placement versus each of the other blast field variables, namely, overpressure, 
dynamic pressure, and overpressure impulse. As expected, we found in each case 
a very iarge point scatter about a trend in which displacement increases with 
increases in any of the blast field variables. We interpret this as follows: 
displacement is causally dependent upon Ig; correlation between displacement 
and the other blast field variables (with very large scatter) is accidental, 
existing only because these variables have a tendency (in the data base) to 
have high values when Iq has a high value (and low values when Ig has a low 
value). 
(U) Upon analysis of the results for the two vehicles and 
two orientations we find that: for fairly small values of Ig turnover occurs 
in the side-on orientation, but not in the face-on, causing measured displace- 
ments to be larger side-on; in the data base there are no large displacements 
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for the M38A1, especially face-on, nor small displacements for the WWII, 
especially side-on; over most of the range of significant overlap of the data 
the disparity between side-on and face-on displacements (for a given value of 
Iq) is about the same for each vehicle. However, the WWII vehicles were ex- 
posed to higher values of Iq than the M38Al's and consequently suffered larger 
displacements. For large values of Iq (and hence displacement) the side-on 
and face-on displacements are close in value. This is because turnover of the 
vehicles (in both orientations) takes place, almost certainly, early during 
the blast pulse. This tends to randomize vehicle orientation, thus making the 
blast force on the vehicle over the remainder of the blast wave and the effec- 
tive resisting force over almost all of the trajectory approximately the same 
for the two (initial) orientations. (We expect that this would be true for 
any initial orientation. ) 


3.1.2 (U) Calculations of Vehicle Motion 


3.1.2.1 (U) Translational and Rotational Velocities 

(U) We made several calculations based on vehicle character- 
istics in order to explain the observations and differences between the 
responses of the two vehicles. Most of the vehicle characteristics we require 
are given in Reference 2. Thése are given in the following table along with 

our estimates of the area, A, presented to the blast by each vehicle (based 
on overall sketches, not detailed drawings of the vehicle), and our calcula- 
tion of the amount of inertia about the overturning axis for the side-on orfen- 
tation, lore (obtained by use of the parallel axis theorem and the listed values 
of the moment, Ioy> about the center of mass). 


VEHICLE/ HEIGHT HALF WIDTH MASS UF I I A Us 
ORIENTA- OF CM OF VEHICLE VEHICLE,M CM OT COEFF. OF 


2 2 
TION (cm) (cm) (gms) (gm cm) (gm cn) (cm?) FRICTION 
64 


WWII, SO 70.2 1.09x10® 2.54x109 1.25x10!9 2.8x10 0.70 

WWII, FO 64 1.09x10® 1.4x104 0.25 | 
M38A1, SO 67.6 71.6 1.19x108 3.00x102  1.46x10'? 3.0x104 0.70 | 
M38Al, FO 67.6 1.19x106 1.5x104 0.25 

(U) The listed characteristics indicate that only small 


differences are to be expected between the responses of the two vehicles. For 
example, the M38A] has a slightly larger mass and overturning moment about the 
center of mass than the WWII: but it also has a slightly jarger area presented 


27 
UNCLASSIFIED 


7 See ee eee Se et Pe ee ee, tee 2 ee ee Le oe?’ So 2 de Se 8 8 EO Se es SF. ES a ee Oe Be: ee ne ne ee ee a se be hs ee 


fore, the differences in the two vehicles tend to cancel one another both for 
translational and rotational response. 


(U) Using the above data we first performed a few very simple 
calculations applicable to both vehicles. 


(1) Taking Iq = 80 kPa-sec = 8 x 10° dyne sec/cm* we find, using 
MV. = IgA where M = vehicle mass and Vp = initial velocity 
Vo 72x 104 cm/sec 2 450 miles/hour. 


(This is for a side-on orientation.) In the experiments subjection of the 
vehicles to blast with Iq considerably less than 80 kPa-sec caused the vehicles 
to be demolished, as would be expected with such a value of the (calculated) 
initial velocity. 


(2) Taking Iq = 10 kPa-sec, we obtain 
Vo ¥ 2.5 x 103 cm/sec * 57 miles/hour. 


If we assume that most of this impulse is delivered in 0.1 sec, the vehicle 


acceleration will be 2.5 x 104 cm/sec* and it will travel a distance of 1.3 
metres. , 
(u) "In general, such calculations show that displacement of 


the vehicle during the intense part of the blast wave is fairly smal] compared 
with the measured displacement of the vehicle. 


(3) The average force acting on the vehicle during delivery of the 
impulse in (2) is: 


5 2 
10" syge secfem” x 3 x 104 om? = 3 x 10'° dynes 


This is much larger than the force of friction, which is 


UNCLASSIFIED - 
to the blast direction and a slightly larger overturning moment arm. There- 
| 
| 


Ferric = 1.1 x 106 gms x 980 dynes/gm x 0.7 


= 7.8 x 10° dynes. 

(4) If we take Ig = 1 kPa-sec delivered in 0.1 sec, the average force 
on the vehicle is 3 x 109 dynes which is still considerably larger than the 
opposing frictional force. Thus even for the smallest blast impulses in the 
experiments the force exerted by the blast wave on the vehicle is greater 
than the opposing frictional force. (For the face-on orientation the vehicle 
are> is on. bout one-half as large but the coefficient of friction is only 


28 
UNCLASSIFIED 


eS ott eat ite. 
CT MOAN es 
i . SRK KRRR i 


_ 2 


UNCLASSIFIED 


slightly over one-third as large as in the side-on case, so the statement js 
true in both vehicle orientations. ) 

(5) We made similar calculations for overturning. Here we assumed that 
the effective moment arm for the force on the vehicle due to the blast is the 
height of the center of mass. Then with Iq = 10 kPa-sec, time of impulse 
delivery 0.1 sec, we find the angular acceleration, a, to be 
2 


a ® 150 radians/sec 

After 0.1 sec, the angular velocity, w, is 
w * 15 rad/sec = 2.4 revolutions/sec 
and the angular displacement, 9, is 
6 = 1/2 at® = 0.76 radian # 43°, 
Here we have ignored the torque owing to gravity, Tg» tending to prevent over- 
‘ 
) 
| 
i] 
' 
| 
: | 
' 
' 
, 
: 

| i 
: 
i 
Lt 
q 
| 


turning. Again using the tabular data, we find 


Tg = 1.1 x 108 gms x 980 dynes/gm x 70 cm 


mW 


7.5 x 1910 dyne cm 
The torque applied by the blast wave is [see (3) above]. 


3 x 10! dynes x 70 cm * 2 x 10!2 dyne cm, 
which is much larger than Tg: 


(U) In general, the angular displacement during the blast 
pulse is much smaller than the angular displacement well after the pulse (as 
measured wherever the vehicle is found overturned) just as was found true for 
linear displacemert; when the bulk of the blast impulse is received over a 
substantially larger time period, 0.5 seconds, for example, the linear and 
angular displacements during the pulse time become quite significant. 

(6) From the calculation in (5) we see that for Iq = 10 kPa-sec the 
vehicle certainly overturns (side-on); for Ig = a few kPa-sec, it certainly 
overturns but for Ig ~ 1 kPa-sec it appears marginal and we cannot draw any 
conclusions from such a calculation. This is because as the vehicle starts 
to rotate a larger area is exposed to the blast wave and the effective moment 
arm also increases. Also, the moment arm associated with the force of gravity 
decreases and when the vehicle has rotated approximately 45°, the torque due 
to gravity reverses direction and helps to overturn the vehicle. With Iq = 
] kPa-sec, delivery time 0.1 sec, the angular displacement at 0.1 sec is 
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owing to gravity rapidly diminishing. 

(U) In the experiments, vehicles exposed side-on to a dynamic 
pressure impulse ~ 1 kPa-sec sometimes overturned and sometimes did not. 

(U) The above simple calculations are applicable to both the 
WWIT and the M38Al vehicles, there being very little difference between the 
two. 


UNCLASSIFIED 
4.3° and the angular velocity is 86° per second, with the restoring torque 
N 
; 
t 


3.1.2.2 (U) Models of Vehicle Motion 
(U) Since the Ig versus displacement data fits for the two 
vehicles appear to be somewhat different and only partially accounted for in 
the previous discussion we attempted a further calculation. The forces we 
should expect to act on the vehicle during its motion, other than the blast 
force, are: a frictional force — during the part of the motion in which 
the vehicle is sliding over (or gouging) the earth; an air drag force which 
at least in the later stages of the motion would resemble the usual type of 
air drag — vehicle penetrating and imparting some of its momentum to the air 
in its path. Earlier in the motion when the long tail of the blast wave 
involves air going past the vehicle this would not seem to apply and as tine 
progresses the air density in front of the vehicle will vary considerably as 
the air refills the region partially emptied of air by the-blast. The Iq 
versus displacement plots contain data from a considerable variety of shots; 
the air density in front of the vehicle depends not only on ti: vehicle's 
(instantaneous) position but also on the weapon yield, height of burst, and 
ambient conditions. Notwithstanding these difficulties we attempted to write 
a (composite) "equation of motion" for a vehicle and determine an effective 
resistive force in agreement with the displacement versus Iq curve (for each 
vehicle, in each orientation). 
(U) We assume the vehicle responds truly impulsively to the 
blast wave; that is, the blast wave imparts a velocity to the vehicle but the 
vehicle moves only negligibly during the blast wave. We also assume that 
during this time the force resisting vehicle motion is negligible. While 
somewhat idealized, these assumptions are in general agreement with the simple 
calculations already presented. We further assume that the resisting force 
(after the dynamic pressure impulse has been received) is, on the average, a 
function only of the vehicle's instantaneous velocity, 7.e., we ignore fluc- 
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tuations associated with the vehicle's orientation (reiative to both its direc- 
tion of motion and the ground). This appears reasonable from one standpoint: 
the vehicle's instantaneous velocity is seemingly the only continuous variable 
during the motion to which the effective resistive force can be related. Then 
with M, V, A representing the vehicle's mass, velocity, and initial projected 
area on the blast direction, we have 


dV dV dx V 
Me” Max de Wax * FOV) = = ty 1) 
where the resisting force or equivalently aW) 4 a function to be determined. 
Then 
- Sy g(V)dV = ip s where D = displacement (at end of motion) (2) 
) 
But D= crg from our curve fits* (3) 
and = MVq * IgA (4) 


* Here 2 is dimensionless since 


fn D=mC+2 en Iq 


- d(én D) _ dD/D 
ee eatery * aT, 
q q/*q 
In CGS units, C has dimensions: 
[c] = Com}® *! [gm)~* [sec]® 


To obtain the values of C for the two vehicles/orientations from the fits 
listed in Table 2.4, where the fits are of the form D = AIG with D in metres 
and Iq in kPa-sec, we have 


. 1OOA 
(104)4 


: 4 dynes sec 
since Jm = 100cm and 1 kPa-sec = 10 cymes see 
cm 
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. where Vo is the vchicle's velocity immediately after receipt of impulse Iq. 
M k 
So D= C(x Vo) ; (5) 


Differentiating the integral with respect to its lower limit, Vo, and substi- 
tuting this relation for D on the right hand side yielas: 


“Big Vy SYIAV = aye U5 HI age Uh 10] * ate (6) 
L L L-] 
Thus (Vo) = ay UW COR Yo) J= RR) Vo (7) 
L -] 
So a(v) = & (ty ev (8) 


and the resisting force is given by the right hand side of the equation 


dV Vv 


£ 
. . WAYS Ay 182 
Mat ~gvy Yc Hw) 2 


(9) 
L 2-2 
Finally, me = gl = - 4h) (10) 
(U) Compariser of the resistive force, r.h.s. of Equation (10), 


which contains y2-e with the values of 2 given in Table 2.4; we see that the 

velocity dependence of the effective resistive force varies greatly and that 

2-2, the exponent of V, can be either vositive or negative to accord with the 
data fits. Thus for low values of V, the resistive force approaches zero for 
2-2 positive, infinity for 2-2 negative, neither of which is physically real- 
istic. This sort of occurrence can perhaps be expected with empirical data 


fits; it does not mean that the effective resistive force derived is necessarily 
invalid over substantial portions of the vehicle's motion. With the resistive 
force derived we can calculate velocity and displacement of the vehicle as a 
function of time by integrating the “equation of motion". However, we have no 
physical interpretation for the derived resistive force so we have only used 

it as a partial check on physically motivated models. 

(U) The displacement versus dynamic pressure impulse curve fits 
contain several phenomena. For small enough values of Iq the vehicles will 
slide across the ground opposed by a conventional frictional force. This 
occurred for both vehicles in the face-on orientation and for the M38Al side-on. 
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The WWII in side-on orientation was not exposed to blasts with such small values 
of Iq. For larger values of Iq vehicle turnover occurs quite readily in the 
side-on orientation, less readily face-on. The M38Al was not exposed to any 
blasts with sufficiently high Iq for turnover face-on. As Iq ts raised and 

the vehicle initial velocity increases air drag can become important, at least 
in the later phases of the motion when the vehicle pushes air which has moved 
back in toward ground zero. At sufficiently high values of Iq the vehicle can 
be partly or completely lifted off the ground, so that there is no ground fric- 
tion term. If there is no air drag term either, the vehicle is in free flight, 
prebabiy tumbling in some fashion, but is eventually pulled down and strikes 
the earth. Thus sliding, rolling, bouncing across the ground along with goug- 
ing of the ground all can take place at one or another time during the motion. 
This, along with the variety of weapon yields, burst heights, ground ranges, 
and ambient conditions, is the reason that the best that can be hoped for is 

a "composite equation of motion". 

(U) With these considerations in mind we consider two cases: 
(1) low values of Iq such that only sliding takes place; (2) higher values 
where the other mentioned phenomena also take place in which we include a 
constant frictional resistance along with a conventiona’? type air drag resis- 
tance. There are clearly some portions of the motion for which these assump- 
tions are inadequate. Consider first the case of small dynamic pressure impulse 
so that only ground friction is important. We take the equation of motion to 
be (it is a legitimate equation for this very restricted case): 


: 


4 
N 


MSE =. msl = - where kK is the constant (11) 


frictional resisting force. Integrating this equat‘on yields 


eo ew EE a ee 


D=MVo/2K (12) | 
The coefficient of friction, u, is then 
u = K/Mg (13) 


For the M38Al in side-on orientation, we have for Ig= 1,0 46.6 cm; with 
values of Ig < 1 kPa-sec, the M38Al did not turn over. 

(U) Using the tabulated data of vehicle characteristics (A 
and M values) we then obtain 


u = 0.71 
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Such good agreement with the measured value of 0.7 is probably fortuitous. 

(In fact over the range of the data 0.597 £ Iq < 1.04 kPa-sec where vehicle 
turnover does not play a role, such calculations provide values of y between 
0.68 and 0.83.) 

(U) For higher values of Ig we include air drag in a "compo- 
site equation of motion" in order to check the result against experiment. Our 
starting equation is then: 


Moy = mvs = -(K) + Kev2), which upon integration and a (14) 
little simplification yields 
K 
M 2.9 
D = in(1 +7 Vo) . 15 
2K ( K, vo) (15) 


Since MVo = I,A, Equation (4), we can simply take two values of Iq in the 

desired range, calculate V, and find D from our curve fits, then use Equation 

(15) to determine K; and Ko. However, the results are quite sensitive to the 

choices of I, values. We found, however, that by fitting the data (for a given 
vehicle/orientation) over a limited range to the expression for D (somewhat 

crudely — we have no convenient least squares program for such a functional 

dependence) we can obtain reasonable results. [That this expression for'D- 
cannot fit the data over a large range is clear Since it is incompatible with 

the fit D = AIS; the latter fit corresponds to the resistive force previously 


derived, F(V) = M (Aye ye-e, Ecuation (10). ] 
COM 


(U) We note that for small Vo, j.e., small Iq, according to 


Equation (15), 
m ke 2 4 2 
O* 2K OR Vo * 2K Yo > 


the same result achieved with no air drag. For large Vy, 


ee ee ee | 
2Ko Ky 0° 0 | 
or D becomes proportional to 2n I. So Equation (15) can only be fitted to 
very limited ranges of data. 
(U) We fitted Equation (15) to data over a limited region 
for a few cases. Taking the WWII vehicle, face-on, over the region 1 < Iq < 
8 kPa-sec we determined reasonable values of K, and Ky [such that Equation (15) 
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predicts displacement to within 15%]: 


K, = 2.35 x 108 


Kp = 850 
K 
Th = = Q, 
en u He 0.22 


which is in reasonable agreement with the measured value, u = 0.25. 

(U) Computing K2Vp@ we find that for Ig = 8, K2V@, the air 
drag term in the (composite) model, is almost 4 times as great as the fric- 
tional term (initially). For I, =, KoVoe << Kj so that here the air drag 
term is negligible throughout the motion. 

(U) We also used Equation (15) for the M38A1, face-on, over 
the range 1 < Iq < 4 kPa-sec (about the range of tl data for this case). 
Following the same procedure we obtained a value for the coefficient of 
friction » = 0.30 which is to be compared with the experimental value 0.25. 
In these cases there were no vehicle turnovers. 


(U) Finally we used Equation (15) for the WWII, side-on, over 
the range 1 < Iq < 8. We obtained the result u = 0.15. : 
(U) This 7s much smaller than the measured value of 0.70, 


corresponding to the fact that all of the vehicles turned over, almost all 
180°; thus the vehicles rolled, bounced, and may have partially or completely 
lost ground contact for part of the motior:. Thus the measured coefficient of 
friction is much greater than the physically effective value. 

(U) Again, calculating KoVo* we find that for Ig = 1 the air 
drag force is (initially) a little less than one-half the frictional force 
while for Ig = 8 the force attributed to air drag in this model is initially 
27 times as large as the frictional force — clearly an important effective 
force over a significant portion of the vehicle motion, even if not a conver- 
tional air drag force. 

(U) Finally, each of the "composite equations of motion", 
Equations (10), (11), and (14), can be integrated enabling us to obtain 
vehicle velocity and displacement as functions of time. (However, we have no 
experimental data with which to compare such model-predicted values.) In 
general, we find that the vehicle motion takes place over a time period ~ 

one to a few seconds. Thus in most cases, al] three models predict that 
vehicle motion persists considerably longer than the intense portion of the 
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blast wave, but in some cases not longer, or not much longer, than the posi- 
tive phase duration. Therefore, in some cases the air drag in our model 
corresponds at least roughly to a conventional air drag term (with varying 
air density) while in other cases it may bear little resemblance to a con- 
ventional air drag force. 

(u) In sum, the various calculations described provide some 
insight into the meaning of the displacement versus Ig curves. Further, they 
indicate that differences in results for the two vehicles are not explained 
to any significant degree by differences in the physical characteristics of 
the WWII and M38Al vehicles. 


3.1.2.3 (U) Surface Effects and Explanation of Remaining Discrepancies 
Between Results for the WWII and M38A1l Vehicles 


(U) Most of the differences in results for the two 1/4 ton 
trucks have already been explained, the causes for the differences being: 
exposure to different ranges of I, values; relatively small regions of over- 
lap of the data for the two vehicles compared with the total data range for 
each vehicle; fairly large scatter of the data points, especially for the WWII 
vehicle, this in turn being caused by the variety of shots which give rise to | 
the data for this vehicle — there being a small, but not completely negligible 
height of burst effect present in the data. We note that the curve fits, good 
straight lines on the 1eg-log plots over the range of data for each yehicle/ 
orientation, would not necessarily be well represented by such straight lines 
if the data extended over a larger range — as it does when we combine the 

data for the two vehicles, which it now seems reasonable to do since differ- 
ences in the data for the two vehicles appear to result from causes other 

than differences in the physical characteristics of the two vehicles. 

* (U) The previously described calculations and data compari- 
sons Ted us to look carefully at the displacement versus Iq data for tne two 
vehicles in the region of overlap of the respective data ranges (in each 
orientation) and to associate the vehicle overturn data with the results 
while making the comparisons. For the two vehicles in face-on orientation 
no vehicle overturns occur in the region of overlap of the data. In this 
region there are only a few data points for each vehicle, 10 for the WWII 
and 7 for the M38A1. There is considerable scatter. The combination of 
few data points and large scatter causes considerable difference in least 
squares fits to the data for the two vehicles in their common region. 
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However, the data for the two vehicles appear to be quite compatible in this 

region; there certainly is no large difference. (The overlap region only 

spans 1.3< Ig <4 kPa-sec. ) 

(U) In side-on orientation there are many more data points 

for each of the vehicles in the region of overlap of the data. The range of 

overlap is 1 £1, < 8 kPa-sec. The slopes of curves for the two vehicles in 

this region are quite different. However, the only significant difference in 

the data sets lies in the region 1 < I, < 3 kPa-sec. Now as we have seen I 

~ ] kPa-sec is marginél for vehicle turnover. The WWII vehicle was not exposed 

to any blasts with Iq <1. The M38A1 data contain 5 points with I, < 1.04; 

in 4 of these cases the vehicle did not turn over. The 5th instance, Iq = 

1.04, D = 0.3\1 appears peculiar in that a 90° turnover occurred — and a 90° 

turnover results in a measured displacement of about 1.4m. In most cases 

examination of our data shows, for small and moderate displacements, that the 

displacement associated with turnover (90° or 180°) accounts for a large 

part of the measured displacement, the latter almost always being somewhat 

larger than the turnover-caused fraction of the measured displacement. (This 

statement applies to side-on orientation; there were no turnovers face-on with 

only small or moderate displacement; when a very large displacement occurs, in , 
either orientation, the vehicle may turn over 360° or more_which is not 

evident in our data — or it may be dexiolished, which is evidert.) 

(U) In tracing back the difference in responses of the two 

vehicles in the region 1 s Iq < 3 kPa-sec, we found first that the M38Al | 
turned over 90° whereas when the WWII turned over it was turned 180°. There were a | 
couple of instances with Iq ~ 3 kPa-sec in which the WWII vehicle did not 

turn over; each occurred with a low yield shot so that, because of the diffrac- | 
tion effect previously described, the displacement was considerably larger | 
than occurs for a similar Iq value associated with a higher yield shot (longer 
blast wave). A 180° overturn causes twice the displacement as a 90° overturr 
and, since most of the measured displacement is caused by overturning, the 
WNIT exhibited larger displacements in the region 1 < Ig < 3 kPa-sec than the 

M38A1 in all cases (owing either to 180° overturn, as opposed to 90° for the 

» M38A1, or to a low yield diffraction effect). Upon further inquiry we found 


; that the WWII vehicle had been exposed to the blast while standing on loose 
' soil while the M38Al stood on firm ground. We believe this accounts for the 
180° WWII overturns compared with the 90° M38Al overturns for the following 
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reason. As the vehicle starts to slide (and rotate simultaneously) the loose 
dirt will be pushed ahead of the tires. When the vehicle has rotated 90° the 
top of the side of the vehicle farthest from ground zero strikes loose dirt 
pushina it in somewhat. Thus the vehicle actually rotates somewhat more than 
90° in a single overturn, the tires elevated a little witn a pile of displaced 
loose dirt under them and the top of the vehicle side farthest from ground 
zero depressed a little. This makes it easier for the vehicle to overturn 
again than if it were lying horizontal. Therefore, over a small range of Iq 
just above Iq ~ 1] kPa-sec, we should expect the WWII, placed on loose dirt, to 
be likely to overturn 180° while the M38A1, placed on firm soil, is likely to 
overturn only 90°. 

(U) This explanation, qualitatively, appears capable of 
accounting for the only as yet unexplained difference in the data for the 

two 1/4 ton trucks. We also note in this connection that the results for the 
2-1/2 ton trucks, Figures 2.5 and 2.6, and Table 2.4 ar? very similar to those 
for the M38Al 1/4 ton trucks. The 2-1/2 ton trucks also stood on firm soil. 
(U) As we noted in a previous section, the displacement 
versus Ig plots and fits are based on gage-derived Iq values for both near- 
ideal and non-ideal shots. As discussed elsewhere in this report we have been 
able to improve upon the Iq estimates for non-ideal shots by use of displace- 
ment data (the vehicle in essence playing the role of a gage, after being 
properly calibrated). If we were to use the best Iq estimates we have obtained 
for plots of displacement versus Iq we would find not only a rather modest 
improvement but also that the discrepancy between data for the two vehicles 

is somewhat diminished where the data scatter is large. This tends to 
corroborate the finding that there is no significant difference between the 
two vehicles and that our explanation of the differences in the response data 
is correct. 

(U) Finally, since we have found that the differences in 

data for the two vehicles are due to factors other than differences in the 
physical characteristics of the two vehicles, it is meaningful to combine 

the data for the two 1/4 ton trucks, in each orientation, and to obtain curve 
fits for the combined data. 


3.1.2.4 (U) Final Dynamic Pressure Impulse versus Displacement 
Least Squares Fits 


(U) We combined data for the two 1/4 ton trucks and obtained 
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fits for Iq versus displacement and displacement versus Ig as shown in 
Figures 3.1 and 3.2. When the cata are combined the ranges of displacement 
and Iq (in each orientation) are greater than for either vehicle alone and a 
straight line log-log fit is not as good as a quadratic; there is a signifi- 
cant curvature. The statistical measures associated with the fits are given 
in Table 3.1. 

(U) We likewise combined the data for the 2-1/2 ton trucks, 
REO M35 and GMC M135. There is insufficient data for a side-on fit for the 
GMC M135 and no face-on data for this vehicle. Therefore, we obtained data 
fits for the REO M35, side-on and face-on, and a combined fit to the side-on 
for the two vehicles. (The few data points for the GMC M135 do appear to 
have slightly higher I. values for a given displacement than is the case for 
the REO M35.) 

(U) We obtained both linear (log-log) and quadratic fits to 
the combined data; the latter is the better and is shown in Figure 3.3. The 
results are shown in Table 3.1. The results for REO M35 alone, side-on and 
face-on, are given in Figures 2.5 and 2.6 and Table 2.4. 

(U) For the combined displacement-Ig fits for the two 1/4 
ton trucks, and for all of the displacement-Iq fits for the 2-1/2 ton trucks,’ 
we used the iterated (improved) Ig values previously discussed for the fon- 
ideal shots. 


q 


3.2 (U) DISPLACEMENT VERSUS SCALED GROUND RANGE. Results are given in 
Figures 3.4 through 3.11 and Table 3.2. 


3.2.1 (U) Scaling 


(U) Scaling ground ranges used here and elsewhere in this report are 
given by 
x = SqR (Modified Sachs scaling) (16a) 
or by xX = SqR ; (16b) 
P, 1/3 
- () 1] 
Here Sq * (7q-7) MEE) (17a) 
P. 1/3 
= (2 a 
and S47 (74-7) yh (17b) 
where Po = ambient pressure, 


a number used in the assumed scaling law, and 
weapon yield (KT). 


n 
W 
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(U) We tried several values of n (in addition to 1/3) for each subcase: 

n = 0.37, 0.4, and 0.43, and in some instances, additional values of n. In most 

of the subcases in this section the optimum value of n was found to be near 0.4 : 
(U) In general, the data have a rather large scatter (since there are 

many uncontrolled variables involved). Thus fairly large errors can be expec- 

ted in inferring displacement from scaled ground range. (See Table 3.2.) 

(U) We also attempted to determine "scaled displacement versus scaled 

ground range" (analagous to "scaled impulse vs. scaled ground range", Section 

3.3). To do this we proceeded as follows. For the scaled impulse case we 


have : 
Silg = f(SgR) (18) 
a, « plo t*) * 273" (4.2, me os a 
where Ss; = Ss 
1 288 Pa yess 


\foairied Sachs Scaling), 


Ty = ambient temperature 
and f indicates a functional relationship. 


; aire | 
Then Ig * gy f(SaR) (20) 


- 


and since displacement, D = ale (see Section 3.1.2.2 and/or Appendix C) we 
have Ds Als> f(SgR)] . This can be rewritten 
stp = aLF(syR)I” , (21) 


where the left hand side now represents scaled displacement. In individual cases, 
ji.e., specified vehicle and orientation, we know A and 2 from our data fits of 
"Displacement versus Dynamiz Pressure Impulse". From out data fits of "Scaled 
Impulse versus Scaled Ground Range" we also have a functional form f containing 
coefficients ag, a], and ap whose values we have determined. If we use these 
results, scaled displacement may be regarded as known and no fitting of data is 
required. However, there is considerable scatter in the data determining A 

and £ (i.e., the data for "Displacement versus Dynamic Pressure Impulse") and 
there is also scatter in the data for determining F(SqR) (i.e., the data for 
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"Scaled Dynamic Impulse versus Scaled Ground Range"). Therefore, we simply 

took sfp = G(SgR) where G(SgR) is a function to be determined by the data and 

it not required that G(SgR) = ALf(SqR) 1° (the right hand side being already 
determined by the above mentioned data fits). This expresses "scaled displace- 
ment" as a function of scaled ground range. 

(U) We fitted scaled displacement, stp, to the functional form of 
Equation (2), Appendix C, with x; = x = scaled ground range and x9 = x2 (in two 
cases where we have sufficient data for such a fit: WWII Vehicle, Ideal Sur- 
face Burst with Yield > 1 KT, side-on and face-on vehicle orientations). We 
regarded 2 as fixed. The resulting fits were not quite as good as the fits 
obtained without scaling displacements. It is, of course, possible that by 
also varying £ the results could be improved. This, however, appears unwarranted: 
if the value of 2 to be used in an application of these results cannot be that 
obtained from "Displacement versus Scaled Dynamic Pressure Impulse" results, 
scaling the displacement has no utility. Therefore, we pursued this no further. 


3.2.2 (U) Effect of Dust 
(U) Having obtained plots and data fits of displacement versus scaled 
ground range for various burst conditions we compared the results. (The group- 
ing of the shots according to burst conditions will be discussed in the next 
section.) Comparing Figure 3.5 with 3.4 (surface burst) we see that the effect 
of dust over a limited span of scaled ground ranges is to greatly enhance the 
displacement. (Since n = 0.4 for both of these figures scaled ground ranges 
are directly comparable. Both figures pertain to the WWII vehicle side-on.) 
At very short range displacement is enhanced by a factor of 6. At long range 
the dust-laden blast wave is weaker than the ideal since energy is sapped from 
the wave by the dust and the dust ultimately falls behind the wave, no longer | 
contributing significantly to the dynamic pressure impulse. 
‘(U) Similarly comparing Figure 3.7a with 3.6 (air burst) we again see 
that the effect of dust is to increase displacement by up to a factor of 6 : 
over a limited spatial region. In this case n = 0.43 for both figures and 
the vehicle involved is the M38Al in side-on orientation. 
(U) Finally, comparing Figure 3.11 (light dust) and 3.10a (near-ideal 
air burst) we obtain the same result for dust enhancement of displacement, 
this time for 2-1/2 ton trucks (with n = 0.43 for both figures). 
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3.3 (U) DYNAMIC PRESSURE IMPULSE AND SCALED DYNAMIC PRESSURE IMPULSE 
VERSUS SCALED GROUND RANGE 
(U) Results are given in Figures 3.12 through 3.21 and Tables 3.3 and 3.4. 


3.3.1 (U) Surface Burst 
(U) We first discuss the surface burst data, ideal/near-ideal blast 
data spanning the yield range of 20 tons - 3.5 MT. Here, as in the case of 
displacement versus scaled ground range, we consider both scaled and unscaled 
dynamic pressure impulse and use both 


x= SqR and x= SR . (16a,b) 


In the latter case we again took values of n (in addition to 1/3) of 0.37, 0.4, 
0.43, and in some cases additional values — somewhat higher values when data 
over a very large range of yields were fitted. 
(U) For both scaled and unscaled dynamic pressure impulse we obtained 
two sets of fits to the experimental data, one using only nuclear shots ranging 
in yield from 39.5 KT to 3500 KT and one including these shots along with HE 
shots extending &:-: ,: |< range downward to 20 tons. We consider first of al] 
the nuclear case. 

d.9.1.1 (cv) Nuclear Events : 
(U) For both scaled and unscaled impulse we have two fits, one 
extending over the entire scaled ground range and one over ‘the large scaled 
ground range domain of the data. See Figures 3.12a, 3.12b, and 3.13. In 
each case the data fit over the entire range reaches a minimum near the end 
of the data range (and would thus be incapable of extrapolation). Comparing 
the fits over the entire data range we see that (Tables 3.3 and 3.4) the statis- 
tical measures of goodness of fit are a little better for scaled than for 
unscaled impulse even with optimum value of n, 0.425 in this case. [The 
ground range scaling factor, Sy depends upon n. See Cquation (17b).] The 
results with n = 0.4, however, are not greatly inferior to those with n = 0.425. 
Comparing the statistical measures for the scaled and unscaled impulse fits 
we have: scaled, Srn 1,C * 0.1637, oc = 0.1611 and Ey = 15.11%; unscaled 
corresponding values 0.1756, 0.1774, and 16.64%, respectively. We note that 
this is in spite of the fact that 1-R2 is 0.0120 for the unscaled case and 
0.0173 for the scaled case. Thus R@ cannot be used to determine which is the 
better of two fits when one or both variables are different in one fit than in 
the other. 
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(U) Comparing the two fits over the large scaled ground range 
end of the data range, again with n = 0.425 determining Sa we see that in 
this case the statistical measures of goodness of fit are a little better for 
unscaled than for scaled impulse. (See Tables 3.3 and 3.4.) In each case 
the minimum of the fitted function here is well beyond the range of the data 
so that reasonable extrapolation is possible. 
(U) In obtaining data fits over the large x end of the data 
region there is some arbitrariness concerning just how much of the data range 
to include. There is in each case substantially more data in the low and 
moderate x than in the large x portion of the data range. In each case we 
obtained additional fits covering a greater or smaller range than indicated 
in Figures 3.12a, 3.12b, and 3.13. We selected what we regarded as the best 
fit on the basis of accuracy of fit to the data points used. number of data 
points used and a subjective criterion: that when the fit obtained is coupled 
with the portion of the total data fit covering the remainder (small x) part 
of the scaled ground range domain — the result is a curve over the entire 
range essentially coinciding with a good eye-drawn curve. 
(U) In further analyzing these results we can say that the 
prediction capability afforded is quite good. There is much less data scatter 
than in corresponding displacement plots. This is partly due to the fact that 
more uncontrolled variables affect displacement than dynamic pressure impulse, 
but is also partly due to the fact that in the data reduction process (by 
many persons over a long time period) there has inevitably been some smooth- 
ing of data for the blast field variables as a function of ground range for 
each shot. Also, some of the latter results are based, at least in part, on 
calculations. (This is not the case for displacement values. ) 
(U) Finally, we remark that in neither the scaled nor unscaled 
impulse plots is the scaling (only ground range is scaled in the latter case) 
perfect; i.e., the data for the 4 shots do not coalesce perfectly. In both 
cases the Nectar-6 points lie a little above the Zuni-3 (and curve slopes 
based on each set separately would differ slightly). These points make up 
the upper end of the curves. Similarly Tower down on the curve the Koon-3 
points lie a little above the Lacrosse-1 points in their (small) region of 
overlap. 
(U) Overall there is not much difference between the relia- 
bility of the two sets of fits. Both are quite good; the scaled or unscaled 
impulse results may be used with equal confidence. 
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3.3.1.2 (U) Nuclear Plus HE Events 
(U) We now consider the corresponding sets of scaled and 
unscaled impulse data curve fits over the entire yield range of 20 tons to 
3500 KT. For the unscaled impulse data we found that for the scaled ground 
range x = SqRs the optimum value of n (upon which S4 depends) is about 0.46, 
the larger value of n in this case being required to coalesce the low yield 
(non-nuclear) data on the same curve as the higher yield shot data — causing 
somewhat more separation of the Nectar-6 and Zuni-3 points at the upper end of 
the curve. See Figures 3.14 and 3.15. 
(U) The comparisons for the curve fits for scaled and unscaled 
impulse are very similar to those just discussed over the smaller yield range: 
the statistical measures are a little better for the scaled impulse fits for 
the whole data range fits but a little better for the unscaled impulse curve 
fits over the large x end of the data. (See Tables 3.3 and 3.4.) Again the 
whole data range fits have a minimum near the end of the data range; the 
large x region fits have minima well beyond the data range and are capable of 
reasonable extrapolation beyond the data range. Once again, the data for the 
several shots are not completely coalesced under the scalings used. In addi- 
tion to a moderate separation of the Nectar-6 and Zuni-3 points at the upper 
end of the curves, the Dice Throw (1 KT, HE) points are somewhat above the 
curve except at large x. (The scaled impulse curve contains one more Dice 
Throw point than the unscaled total data fit. As can be seen from Figure 3.14, 
the Dice Throw point at x = 574.7 cannot reasonably be included in the curve 
fit over the total data range without considerable worsening of the fit. 
Thus the total data fit does not extend to this scaled ground range, though 
the fit to the large scaled ground range data includes this point. We have 
no reason to suppose this to be a bad data point; the functional form 
selected for fitting the data, however, does not provide as accurate and 
reliable a predictive capability when this point is included in the fit over 
the entire data range.) At the large x end of the data the Koon-3 points are 
somewhat higher than the other data points. Deviations from the curves are 
modest, however, and considering the large yield range encompassed by the 


né data, the scalings used coalesce the data quite well. Again, good eye-drawn 
u curves would coincide rather well with the combined curves formed by the 
‘ large scaled ground region fits and the non-overlapping regions of the total 
y data curve fits. (At the low ground range end where the curves are very 
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steep, eye-drawn curves could improve on the fitted curves somewhat more here 

than in the all-nuclear shot case previously discussed where the yield range 

covered is much smaller. ) 

(U) Finally, the predictive accuracy is somewhat diminished 

as compared with the all-nuclear, smaller yield range data fits. (See Table 

3.3 and 3.4.) For example, the RMS errors, £, for the scaled impulse versus 

scaled ground range fits over the whole data range and over the yield ranges 

(20 tons - 3500 KT) and (39.5 KT to 3500 KT) are 19.24% and 15.11%, respec- 

tively. The corresponding results for (unscaled) impulse versus scaled ground 

range are: 20.22% and 16.64%, respectively. Similar results hold for the 

other statistical measures. 

(U) Once more from an overall viewpoint there is little 

difference in reliability between the two sets of curve fits to the data. 

Both are good. The scaled or unscaled impulse versus scaled ground range 

curve fits may be used with equal confidence. 

3.3.1.3 (U) Comparison of Scaled and Unscaled Dynamic Pressure 
Impulse Results 

(U) - To show that this is true we compared values of Iq 

obtained by one of the scaled and unscaled fits for a matrix of values of 

yield and ground range. (We ignored effects of ambient canditions here 

which are small and have a negligible effect on the comparison.) We chose 

weapon yields W, of 0.1, 10, and 1000 KT and values of R/W!/> of 150, 200, 

300, 400 and 600 metres. Here R is the ground range and the specification 

of r/wi/3 was simply for convenience in the computation using the scaled Ig 

versus scaled ground range fits. The actual domain of ground range R in the 

matrix of R, W values is 69.62 < R < 6000 metres. We used Figures 3.14 and 

3.15 and the data fits listed thereon. For the selected matrix of values 

we thus calculated 15 values of Iq over the very large span of yields and 

ground ranges using both the scaled and unscaled I, data fits. (The range : 

of R extended into the extrapolated region of the data fits - large x 

regiun.) The average percent difference between the two sets of results 

was 13.6%. ! 
\ 
| 


3.3.2 (U) Non-Surface Bursts (HOB and Surface Effects) 
(U) In section 3.2.2 we compared displacements for near-ideal bursts 
with those for bursts which gave rise to appreciable dust in the blast wave. 
Prior to plotting and fitting this data, however, it was necessary to group 
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the shots. In order to do this we must have some criterion since the amount 
of dust will vary with each shot for each src.ind range. Since a primary 
objective of the present study is to obtain gocd estimates of Tg» especially 
for non-ideal shots, we first plotted scaled Iq versus scaled ground range for 
the varfous shots. We made this choice because the scaled Iq plots involve 
w'/3 an scaled ground range — whereas (unscaled) Iq plots involve W" with n 
variable, at least over a limited range of values. Having done this, however, 
we find the results obtained to be consistent for scaled and unscaled Iq and 
also for displacement versus scaled ground range. The resulting grouping 

tis as follows: "Light Dust" - Grable-10, Moth-2, Turk-4, Bee-6, Apple I-8, 
and Smoky-15; "Heavy Dust" - Met-12 and Apple II-13; "Non-Surface Near-Ideal" - 
Encore-9, Wasp Prime-9, Yuma-4 (and Wasp-1, although we have no Iq data for 
this shot except what can be derived from Iq versus displacement; but al] 
vehicle displacements were very small for Wasp-1, hence inference of Ig is 
very inaccurate). This‘ grouping is somewhat coarse since the data do not 
warrant fine distinctions as to degree of dust. See, however, Section 4. 

(U) This grouping parallels the scaled HOB grouping for the shots 

— see Table 2.1, with one exception: by virtue of its scaled HOB we would 
expect Turk-4 to be in the "Heavy Dust" group with Met-12 and Apple [I-13. 
However, Turk-4 was exploded over a well-vegetated region, thus creating less 
dust than shots over desert at the same selected HOB. As far as dynamic 
pressure impulse is concerned, the dust, aot the scaled HOB, is the relevant 
factor. 

(U) Figures 3.16a, 3.16b, and 3.17 show that the data for the various 
shots, grouped as indicated, fall] into a reasonably compact set. The Bee-6 
points are a little high, perhaps indicating somewhat more dust than other 
shots in the group. These figures are ror "Light Dust". Figures 3.18 and 3.19 
for "Heavy Dust" exhibit somewhat more data scatter. See Section 4. (In the 
case of Figure 3.18 the value of n is almost irrelevant since Met-12 and Apple 

II-13 have scaling factors which are close in value; thus changing n produces 

very little change in the relative positions of the data points for the two 

shots. ) 

(U) Finally, we wish to evaluate the effect of dust on Ig. Comparing 

Figures 3.16a and 3.12b (each of which has n = 0.4) we see that light dust 

increases I_ by a factor of about 1.2 to 1.6 over the region of scaled ground 

range overlap of the data. Likewise comparing Figures 3.18 and 3.12b (n = 0.4 
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for both) we see that heavy dust enhances Iq by a factor varying from about 

1.5 to 3 over the scaled ground range region common to the two figures. 

(U) The effect of dust on I, can also be seen by comparing the scaled 
Iq versus scaled ground range plots. Thus comparing Figures 3.17 and 3.13 

{or 3.15) we see that over the common region of scaled ground range light dust 
increases (scaled) Ig by a factor varying from about 1 to 3. Again for larger 
ground range values (scaled) Iq falls off much more rapidly when dust is present 
than when it is not. Comparing Figures 3.19 and 3.13 (or 3.15) we see that 
heavy dust enhances (scaled) Iq by a factor of about 3 to 5 over a limited 
spatial region; thereafter (scaled) Iq falls off very rapidly with increasing 
ground range. 

(U) Generally comparing Figures 3.21 and 3.13 (or 3.15) we see that 

at large scaled ground ranges Iq falls off more slowly for a near-ideal air 
burst than for an ideal surface burst. For small scaled ground range Iq is 
much greater for a near-ideal surface burst than for a near-ideal air burst; 
however, our data does not extend to very small ground ranges. We note also 
that the scaled HOB's for the shots we have grouped in the near-ideal air 
burst category differ considerably from one another, which probably explains - 
at least partly the lack of good coalescence of the data for the separate shots 
in Figures 3.20a, 3.20b and 3.21. Z 


3.3.2.1 (U) Comparison of Scaled and Unscaled Dynamic Pressure 
Impulse Results 
(U) We used both the scaled and unscaled Iq versus scaled 


ground range plots (separately) in obtaining preliminary scaled Iq HOB contours. 
(In Section 4 we descrihe the method we finally used in obtaining these con- 
tours; however, the effort described here provides a good comparison of the use 
of the scaled and unscaled dynamic pressure impulse results.) Except fur the 
surface bursts each shot provided data points on *ne contours at it: scaled HOB, 
the latter being different for each non-suriace shot. In using the figures 
just discussed we followed the fitted curve for the shot group in general, but 
allowed for deviations of the individual shots from the group curve (since the 
shots do not coalesce perfectly on the group curve). That is, we assume that 
there are real physical deviations of individual shots from the group curves 
(certainly there are various degrees of dust) — the deviations are nct due 
only to incccuracy of our knowledge of Ig. In any given case we cannot know 

to what extent our inaccurate knowledge of Ig and to what extent real physical 
differences account for deviations from the group curve. However, the group 
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be sufficient data to determine a reliable curve for a single shot. The re- 
sults obtained in this way, using separately the scaled and unscaled I, versus 
scaled ground range fits, for obtaining values of SyIq for the HOB contours 
were very close — nearly always within about 2% of one another. 


3.4 (U) BEST VALUE OF n 

(U) When we started fitting displacement and I, versus scaied ground range, 
the first several cases we tried seemed to have their optimum values of n in 
the vicinity of 0.4. As we have seen, comparisons of figures can be made 
directly only when the values of n are the same. So we adopted n = 0.4 as a 
"standard" for comparison purposes. For our scaled Iq versus scaled ground 
range fits, however, the value of n is fixed at n = 1/3 {modified Sachs scal- 
ing). There is a theoretical basis for this, although effect of dust is not 
included; we have seen, however, that it is not useful to scale displacement. 
Therefore, having improved our estimates of I, end obtained fits for a con- 
siderable variety of cases we wish to establish our "optimum" overall value of 
n — one which could be used for cases not treated herein with the expectation 
that reasonably accurate results would be achieved. 

(U) In Tables 3.2 and 3.3 we list the results for many cases we have fitted. 
The values of n are optimum within approximately + 0.01 when only one 
value of n is listed. (In a few cases, as mentioned previously, changing n 
dees not shitt the points very much relative to one another so that the re- 
sults are almost independent of n. This can happen when cnly two or three 
shots supply data for a curve to be ficted*; generally when there are three 
or more snots, however, the results vary quite significantly with n.) In 
several cases we list the results for the optimum value of n and also for n 

x 9.4. Comparison then shows hcw much better the optimum is than that ob- 
tained wit! n = 0.4. A graphical comparison is obtained by comparing Figures 
3,7a, 3.503 3.6a, 3.8b; 3.12a, 3.12b; 3. 16a, 3.16b**; 3.29a, 3.2Cb. The two 


UNCLASSIFIED 
curve served as an excellent guide in al] cases; in many cases there would not 


TW) “¥ When only two shots contribute data it msy also happen that their data 
ranges don‘t overlap in I, or scaled I (or overlap is based on one or two 
data points) so that the data can be “dgalesced" by a single curve within a 
wide range of n (including n = 1/3} 
(U)** The fits in Figures 3.16a and 3. 16b. iviustrate the case in which one statis- 
ticul meajure of goodness of fit is slightly better for n = 0.4, the other 
two slightly better for n = 0.37. See Table 3.3. The optimum value of n 
in this case is probably between thes2 tvn values. 
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figures in each pair differ only in their values of n. Inspection of the 
entire set of results, Tables 3.2 and 3.3, shows that the optimum-value of u, 
overall, is abcut n = 0.42 or 0.43. Most of the cases for which the data are 
best and there are the most data pcints cluster about such a value. In fact 
all of the "good" cases and almost all of the "fair" and “not so govd" ca: es 
Tie within the renge 0.40 < n < 0.46. Perhaps the two best cases, 3.12a and 
3.14, all nuctear near-ideal surface burst and nuclear + HE near-ideal su'face 
burst Iq versus scaled ground range, have optimum values of n of 0.425 avd 
0.46, respectively, the latter comparatively high value of n being needeu to 
coalesce the data over the very large yield range of 0.04 < Y < 3500 KT. As 
would be expected with many data fits, a goodly fraction of which have rather 
few data points with considerable scatter, values of optimum n will vary some- 
what more than is physically realistic. Figures 3.9, 3.10b, 3.20a and 3.20b 
may well be examples of this. 

(U) So while we have no criterion which allows us to determine "the" optimum 
overall value of n, a value of n = 0.42 or n = 0.43 appears to be a very good 
choice and one which can be used with some confidence for a broader cliss of 
data than encompassed herein. References 4 and 9 show similar ‘alues of n 

for blast damage to drag-type targets. 


3.5 (U) EFFECTS OF NON-DESERT SURFACE (WATER, ASPHALT, DESERT ROCK) 

(U) In shot Met-12 gages and vehicles were placed not only on desert surface 
but also on water (over desert}, asphalt and desert rock; on shot Bee-6 
similar data were taken over asplialt. We did not use any of this data in | 
obtaining the fitted curves. We now wish to compare this data with the 
corresponding desert surface results. We start with Iq {and scaled Iq) 
versus scaled ground range. 


3.5.1 (U) Effect on Dynas ‘c Pressure Impulse 
(U) Placing the data pc +s for a water surface on Figure 3.18 and 
3.19 we see that the Iq values tor water are substantially lower than the 
desert values (except for the lowest point). See Figures 3.18-1 and 3.19-1. 
The Iq values for asphalt are lower than those for water over the whole range 


and much lower than desert, a factor of 2 or more. Placing the asphalt data 
for shot Bee-6 on Figure 3.16a again shows the asphalt Ig values to be lower 
than those for desert (3 of the 4 data points lower, fourth about the same). 
See Figure 3.16a-1. On the average the asphalt Iq values are Tower by a factor 
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of about 1.4. Use of Figure 3.17 with the asphalt data produces the same 
result; see Figure 3.17-1. 


3.5.2 (U) Effect on Displacement 
(U) We next compare displacements for the desert and non-desert data. 
Referring to Figure 3.5-1 (side-on) we see that displacements for vehicles on 
a water surface are about the same as for desert (although the Ig values for 
water are somewhat lower than for desert). Use of Figure 3.9-1 provides the 
same result; in this case we compare face-on displacements. Using the same 
figures we see that vehicle displacement on asphalt is much smaller than on 
desert surface. This is consistent with the considerably lower Iq values 
found over asphalt. 
(U) We also find that vehicles placed on Desert Rock are displaced 
about the same as on asphalt — if we assume no difference between the 
M38A1 and WWII vehicles. This is a good assumption since the displacements 
involved here are large (see Section 3.1.2.3). 
(U) In the case of the 2-1/2 ton trucks placing the Met-12 (heavy 
dust shot over the desert portion of the terrain) displacements (Desert Rock) 
on the light aust plot for 2-1/2 ton trucks, Figure 3.11, the most comparable 
data we have — we see that the displacements over Desert Rock are uniformly 
substantially lower than over desert in agreement with the-above result for 
1/4 ton trucks. See Figure 3.11-1. 


3.5.3 (U) Effect on Displacement versus Dynamic Pressure Impulse 
(U) We next compare displacement versus Ig for desert and non-desert 
cases. Comparing vehicle displacements over water with those over desert, 
Figures 2.1-1 and 2.3-1, we see that the differences are rather smail. There 
is not sufficient data to determine whether or not this is significant. 
(U) Comparing displacements over asphalt with those over desert, aguin 
using Figures 2.1-1 (side-on) and 2.3-1 (face-on), we see that there is no 
difference in the face-on orientation, while the only significant difference 
in the side-en orientation occurs for one Met-12 point with Ig value in «he 
vicinity of ] kPa-sec, i.e., the marginal region for turnover. 
(U) Thus as far as displacement is concerned for a specified value of 
Iq, there is not much difference among water, asphalt and desert. (See, how- 
ever, Sectior 3.1.2.3 for effect of surface on vehicle turnover in cases where 
I 


q 


value is marginal for vehicle turnover. ) 
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3.5.4 (U) Effect on Damage 


(U) Finally we come to a comparison of damage between desert and 
non-desert surfaces. The damage data have inherently large scatter. Refer- 
ring to Figure 3.3-1, we see that the results for water differ little from 
those for desert, though one point appears low — this is in face-on orienta- 
tion. Using Figure 3.1-1, side-on, we see that two of the water points lie 
far below the curve (while the other 3 are on the curve in the damage = 10 
region). For one of the two points, the vehicle apparently was supported by 
a dike retaining the water; in any case it did not turn over despite a large 
displacement (12 metres). In the second case, 78 metre displacement, the 
vehicle probably turned over 360° — but was not damaged as it would have 
been had it been bouncing and/or scraping the earth. 


(U) For asphalt, comparing the side-on data with that of desert, 
Figure 3.1-1, we find that 8 of the 9 points are within lo of the curve 
(mostly above); the ninth is a damage = 10 point 30 above the curve, yet 

not very far (to the low displacement side) from other (desert) damage = 10 
points. In the face-on case comparison, Figure 3.3-1, 6 of the 9 points are 
within lo of the curve again mostly above. ‘Two Of the remaining 3 are damage 
= 10 points, one of these not far from other such points. The otner was 
evaluated as damage = 10 owing to fire rather than blast damage; the blast 
damage was moderate. The remaining point corresponds to a vehicle for which 
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part of the damage was caused by debris rather than blast. Overall, then, 
damage for a given displacement appears to average somewhat higher on an 
asphalt surface than on desert surface. (As we have seen, for a specified 
displacement, I_ is about the same on desert and asphalt surfaces, so that 
any dependence of damage on I. over and above that owing to displacement — 
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(U) * The damage data and results are discussed in a separate report, Reference 
17; figures from that report used here follow Figure 3.21 in this report. 
Figures 2.1, 2.2, 2.3, and 2.8 of Reference 17 correspond to Figures 3.1-1, 


e 


Re 


a 3.2-1, 3.3-1, and 3.8-1, respectively, of the present report. The curve 
Ty fit shown is, in each case in both reports, the damage-displacement least 
re squares data fit for vehicles on a desert surface. However, in Reference 
as 17 the (desert) data points used in the fit are shown whereas in the present 
i report the data points shown are non-desert data points (and ave not used 

in obtaining the fitted curve). 
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which is in turn caused by the dynamic pressure impulse, I, — is nota 
factor here.) 


(U) Lastly, comparing damage over a Desert Rock surface with that 
over desert we have recourse to Figures 3.2-1 and 3.8-1. From Figure 3.2-1, 
we see that there are not sufficient Desert Rock data points on which to 
base a judgement; however, the Desert Rock points do exhibit somewhat more 
damage for a given displacement than the desert points. In the case of 
Figure 3.8-1 (2-1/2 ton trucks) the result is more clear cut: there are 

8 desert rock data points, none in the damage = 10 region and all above both 
curve fits, averaging more than 20 above curve fit 1 and 1.50 above fit 2. 
Thus, it is clear that for a given displacement vehicle damage is greater on 
a Desert Rock surface than on desert. 


(U) Thus the pattern is consistent: for a given I. value, displace- 
ment is, at least approximately independent of the surface on which the 
vehicle is placed (provided, for the side-on orientation case, that I 

above the marginal region for vehicle turnover, ~1 kPa-sec); but for a given 
displacement damage is greater on the harder surfaces, asphalt and Desert 
Rock, than on desert and greater on desert than on water (over desert), j.e., 
a soft, mushy surface. This is quite plausible since bouncing, scraping and 
gouging a hard surface would be expected to do considerable damage to a 
vehicle, whereas a soft, mushy surface has a cushioning effect. 


3.5.5 (U) Plumbbob, Smoky-15 Results 


(U) The data for Plumbbob, Smoky-15 differ moderately from other 
desert data: 

(1) As far as I, and scaled I. versus scaled ground range 
are concerned, the data appear consistent with other shot data. See Figures 
3.16a, 3.16b, and 3.17. 

(2) As far as displacement versus I, is concerned, the 
Smoky-15 data are within the data scatter. Smoky-15 data are included in 
these figures. 

(3) For displacement versus scaled ground range, the Smoky-15 


points are somewhat low, both side-on and face-on; we did not use them in 


the fits. (They are comparable with the asphalt points here.) 
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(4) For damage versus dispiacement* we did not use the 
Smoky-15 data in the fits. On comparing the data with the fits, however, 
we found that the Smoky-15 points averaged somewhat below the desert data 
curve in side-on orientation and above the desert data curve in face-on 
orientation. 


(U) There are only a few points in each instance so that it is 
difficult to draw clear-cut conclusions. However, it appears probable that 
the (modest) discrepancies can be accounted for by combination of the usual 
data scatter and a small inaccuracy in our estimates of dynamic pressure 
impulse. (Our estimate of dynamic pressure impulse may be a little high.) 


(U) * Damage results are reported in Reference 17. 
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TABLE 3.5 (U) Non-Desert Surface Data 
UNCLASSIFIED 


VEHICLE GROUND Iq Displacement 
TYPE RANGE _(M) (kPa-sec) (M) da = SURFACE 


TEAPOT, Bee-6, Side-on 


WHIT 549 €.59 20.1 7 Asphalt 
WHIT 610 4.16 5.7 3 Asphalt 
WHIT 701 1.60 3.6 3 Asphalt 
UWIT 77? 1.02 1.7 3 Asphalt 
TEAPOT, Bee-6, Face-on 
WHIT 549 6.53 13.4 3 Asphalt 
WHIT 610 4.16 1.6 7 Asphalt 
WHIT 701 1.60 1.7 2 Asphalt 
WWIT 777 1.02 0.5 3 Asphalt 
TEAPOT, Apple I-8, Side-on 
M38Al 869 -- 2.7 3 Desert Rock 
TEAPOT, Met-12, Side-on 
WWII 610 28.1 3° 10. Water 
Wet 686 ; 19.6 103 10 Water i 
WWIT 762 13.1 157 10 Water 
WWIT 838 é 7.55 78 5 Water 
WWIT 914 3.94 12 7} Water 
WWIT 610 19.0 68 10 Asphalt 
WWII 686 13.7 59 10 Asphalt 
WNIT 762 8.4 23 10 Asphalt 
WWII 838 5.58 14 6 Asphalt 
WHIT 914 1.72 1.0 ] Asphalt 
TEAPOT, Met-12, Side-on i 
M38A1 610 -- 81.0 10 Desert Rock : 
M38Ai 686 - 54 10 2sert Rock | 
M38A1 762 -- 22 9 Desert Rock 
M38AT 838 - 11.8 6 vVesert Rock 
REQ M35 762 ~- 15.3 8 Desert Rock 
REO M35 762 -- 8.1 9 Desert Rock 
REO M35 762 -- 16.9 7 Desert Rock 
REQ M35 838 -- 10.3 6 Desert Rock 
RED M35 914 -< 5.7 6 Desert Rock 
GMC M135 686 -- 16.7 9 Desert Rock 
GMC M135 762 -- 13.1 8 Desert Rock 
GMC M135 838 -- 6.7 5 Desert Rock 
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TABLE 3.5 (U) Non-Desert Surface Data (Concluded) 
UNCLASSIFIED 


VEHICLE GROUND Iq Displacement 
TYPE RANGE _(M) (kPa-sec) (M) da SURFACE 


TEAPOT, Met-12, Fave-on 


WWII 610 28.1 110 10 Water 
WHIT 686 19.6 91.4 10 Water 
WWIT 762 13.1 88 10 Water 
WWIT 838 7.55 8.8 1 Water 
WWIT 914 3.94 3.3 1 Water 
WWIT 610 19.0 71.3. 9 Asphalt : 
WWIT 686 13.7 41.4 10 Asphalt ° 
_ WII 762 8.4 20 10 Asphalt 
* WWIT 838 5.58 4.0 3 Asphalt 
WWIT 914 1.72 0.5 3 Asphalt 
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‘FIGURE 3.4 (U) Displacement versus Scaled Ground Range - 1/4 Ton 
Trucks, WWII, SO, Near-Ideal Surface Burst, n = 0.4 
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FIGURE 3.5 (U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
WWII, Side-on, Light and Heavy Dust, n = 0.4 
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FIGURE 3.5-1 .(U) Displacement versus Scaled Ground Range ~- 1/4 Ton Trucks, 


WNIT, Side-on, Light and Heavy Dust, n = 0.4, Desert Data 
Curve, Non-Desert Data Points 
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FIGURE 3.6 (U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
M38A1, Side-on (Near-Ideal Air Burst), n = 0.43 
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(U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks 
M38A1, Side-on, Light Dust, n = 0.43 
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FIGURE 3.7b (U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
=; M38A1, Side-on, Light Dust, n = 0.4 
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FIGURE 3.8a (U} Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
WWIT, Face-on, Near-Ideal Surface surst, n = 0.43 
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FIGURE 3.85 (U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
WWII, Face-on, Near-Ideal Surface Burst, n = 0.4 
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FIGURE 3.$ (U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
WWIT, Face-on. Light and Heavy Dust, n = 1/3 
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FIGURE 3.9-1 (U) Displacement versus Scaled Ground Range - 1/4 Ton Trucks, 
WWII, Face-on, Light and Heavy Dust, n = 1/3, Desert Data 
Curve, Non-Desert Data Points 
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FIGURE 3.10a (U) Displacement versus Scaled Ground Range - 2-1/2 Ton Trucks, . 
REO M35 and GMC M135, Side-on, Near-Ideal Air Burst, n = 0.43 
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FIGURE 3.10b (U) Displacement versus Scaled Ground Range - 2-1/2 Ton Trucks, 
REO M35, Side-on, Near-Ideal Air Burst, n = 0.45 
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FIGURE 3.11 (U) Displacement versus Scaled Ground Range - 2-1/2 Ton Trucks, 
REO M35 and GMC M135, Side-on, Light Dust, n = 0.43 
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FIGURE 3.11-1 (U) Displacement versus Scaled Ground Range ~ 2-1/2 Ton Trucks, 
REQ M35 and SMC M135, Side-on, Light Dust, n = 9.43, "esert 
Data Curve, Non-Desert Data Points 
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FIGURE 3.12a (U) Dynamic Pressure Impulse versus Scaled Ground Range - Near-Ideal 


Surface Burst, n = 0.425 
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FIGURE 3.12b (U) Dynamic Pressure Impulse versus Scaled Ground Range - Near-Ideal 
Surface Burst, n = 0.4 
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FIGURE 3.13 (U) Scaled Dynamic Pressure Impulse versus Scaled Ground Range - 
Near-Ideal Surface Burst 
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FIGURE 3.14 (U) Dynamic Pressure Impulse versus Scaled Ground Range - 
Near-Ideal Surface Burst, n = 0.46 
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FIGURE 3.15 (U) Scaled Dynamic Pressure Impulse versus Scaled Ground Range - 
Near-Ideal Surface Burst 
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FIGURE 3.16a (U) Dynamic Pressure Impulse versus Scaled Ground Range - 
Light Dust, n = 0.4 


87 


UNCLASSIFIED 


ee ae es . ; 


oe 


: a é } vie n> He ; 
ete aug es Shas Ree te GS ee ee ah pi RRR 4 ne oa ae Ss Rance ee 


tot rer ttt te aL Rea Renee ene ies eer wendy : 


UNCLASSIFIED. 


en Iq * 5.967 = 0-0198x + 0.00000655 x2 DESERT DATA 
R2 = 0.948 CURVE FIT 


ees _Scaled Ground Range (Metres) 


) 100 


~~ 600° 600 O70 


UNCLASSIFIED 


FIGURE 3.16a-1 (U) Dynamic Pressure Impulse versus Scaled Ground Range - Light 
Dust, n = 0.4, Desert Data Curve, Non-Desert Data Points, 
Bee-6, Asphalt 
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FIGURE 3.16b (U) Dynamic Pressure Impulse versus Scaled Ground Range - 
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FIGURE 3.16b-1 (U) Dynamic Pressure Impulse versus Scaled Ground Range - Light 
Dust, n = 0.37, Desert Data Curve, Non-Desert Data Points, 
Bee-6, Asphalt 


90 
UNCLASSIFIED 


free een mene ees ane Hote eA cy 4. i 
Wisth seperate Parwney: SORE eevee Rene eo nee RMR Meee SS eres ex 


FIT: 2n sity = 3.720 - 0.00663x - 9.0000118 x2 
2 s 0.952 


4__e 
3.0. 

—— eet e Apple 1-8 
— SS eee @ Smoky-15 - 


oo a Po Grable-10 


| 


hartabhinmmtin rr 


Scaied Iq {k®a-sec) 


-eceat 


! 


1o°l. > —-- —-—-~ -- “Sealed Ground Range (Metres) -- ~~ ~- moet . 
area [are], 300 AO 500 B00 700 
UNCLASSIFIED 
FIGURE 3.17 (U) Scaled Dynamic Pressure Impulse versus Scaled Ground Ranae - 
Light Dust 
91 


UNCLASSIFIED 


« 


a Teter be eee e es 


Cee ~ ‘ ees Sekt ao oe . ae se age 1, 
SBA 


Sais Pas se te 


a? 


ene wr oe oy 
5 yt sate <4 s te Routes “y 


areas area : wate 
: i! rd tee 
Ty " Roark Aeana mete ON aR ae marge MAGA? 


a0 AMAR BOT PUREE exo 


UNCLASSIFIED | 


FIT: 2n Fe 


q 3-720’ 0.00663x - 0.000118 x° 
r° = 0.952 


i 


‘ : ‘ 
verrerberss tryin 


ee er. 


19; ————----_ Scaled Ground Range (Metres) 7 
0 100 200 300 400 500 600.7290 
UNCLASSIFIED 


FIGURE 3.17-1 (U) Scaled Dynamic Pressure Impulse versus Scaled Ground Range - 
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FIGURE 3.18 (U) Dynamic Pressure Impulse versus Scaled Ground Range - 
Heavy Dust, n = 9.4 
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FIGURE 3.18-1 (U) Dynamic Pressure Impulse versus Scaled Ground Range - 
Heavy Dust, n = 0.4, Desert Data Curve, Non-Desert Data 
Points 
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FIGURE 3.19-1 (U) Scaled Dynamic Pressure Impulse versus Scaled Ground Range - 
Heavy Dust, Desert Data Curve, Non-Desert Data Points 
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FIGURE 3.20a .(U) Dynamic Pressure Impulse versus Scaled Ground Range - 
Near-Ideal Air Burst, n 2 0.46 
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FIGURE 3.20b (U) Dynamic Pressure Impulse versus Scaled Ground Range - 
Near-Ideai Air Burst, n = 0.50 
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FIGURE 3.21 (U) Scaled Dynamic Pressure Impulse versus Scaled Ground 
Range - Near-Ideal Air Burst 
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(U) SECTION 4 
(U) HEIGHT OF BURST CHARTS FOR I 


q 
4.1 (U) SCALED DYNAMIC PRESSURE IMPULSE VERSUS SCALED GROUND RANGE CURVES 
(U) Prior to constructing the Height of Burst (HOB) Charts, we again 


went through an iterative procedure to obtain Iq values as a function of range 
for each nuclear event (HE information was not used in constructing the HOB 
Charts). The procedure used was similar to that described in Section 2.2.4, 
except that we used the quadratic least squares fits (see Section 3.1.2.4) 

for the combined WWII and M38A1 (1/4 ton trucks) and for the combined M35 and 
M135 (2-1/2 ton trucks) to obtain Toy . The is data for each nuclear event 


ns are shown in Tables 4.1 through 4.3 a derived from Desert Rock or from 
Pas water surfaces were not used in these Tables). Tables 4.1 and 4.2 show the 


data for ideal/near-ideal blast environments for surface bursts (4 events) 

and for air burst (4 events), respectively. Table 4.3 shows the data for 
non-ideal blast waves (8 events). Igg in these Tables represents the original 
estimates (non-iterated) for the ideal/near-ideal events and the first itera- 
tion for the non-ideal events. Tgavg represents the calculated values from 
vehicle displacements (both 1/4 ton and the 2-1/2 ton trucks). - The Iqg -and 
Igavg were then averaged to produce Ig (actual). Note that in some cases 
there were no derived Toavg values. In these cases, vehicles were not loca- 
ted at these ranges so only Iqg was used and these were then labeled as Tg 
The Tq along with range were scaled to 1] KT sea level conditions using modi- 
fied Sachs scaling. The scaled Ig versus scaled range data for each event 
were grouped according to their scaled height of burst. Our procedures 
described in this Section on grouping and applying least squares fits to the 
data for various burst conditions are similar but a refinement to that 
described in Section 3.3. 

(U) The ideal/near-ideal surface bursts were comprised of four 
events: Koon-3, Nectar-6, Lacrosse-1, and Zuni-3. Least squares fits were 
made to these four events. Of the four ideal/near-ideal air bursts, there 
is some uncertainty about Yuma-4, i.e., whether it was ideal/near-ideal or 
non-ideal. We previously considered this to be an ideal/neav-ideal blast 
event since the waveforms displayed steep-rising shock fronts (one of the 
measures for identifying blast wave categories). However, in our final 
resolution of I, for Yuma-4 these values are all above our initial estimates 


q 
105 
UNCLASSIFIED 
ate at a ata 4 oe TS betes SOLA TAY ee antes 
SARTRE ae setbeiatcht onesaeeee ies uate 


DUR ee ee ere Son gE aN Gy bag Ie See PA a POR, ey ge ee eho Se ee ec ee ee ee ee ee ee 
‘i 

4 

Y 


UNCLASSIFIED 


(see Appendix A), implying that these higher values of I, could ve the result 

of precursor action. But Yuma-4 was a fractional KT device and iiffraction 

could contribute to the motion. The derived Tgayg from vehicle Jisplacements 

would then indicate higher values. Data on Yuma-4 is insufficien* to be able 

to distinguish the discrepancy between diffraction effects and precursor 

action. Another questionable event is Wasp-1. Only peak pressure values 

were obtained on this event so the Tao values for Wasp-1 were obtained by 

scaling Encore-9 to the actual Wasp-1 conditions. The scaling using modifie« 

Sachs is straightforward and there is no reason to doubt the validity of the 

Igy values on Wasp-1. However, there are two factors which contribute to the 
questionable data on Wasp-1. First, the actual ground zero was somewhat dis- 

tant from the intended ground zero causing the orientation of the vehicles to 

be at some angle other than side-on or face-on. Our Iq versus displacement 
correlations are either for side-on or face-on orientations. The second 

factor was that the displacements were fractions of a metre. Both factors 

would contribute to large errors in computing laavg from displacements. 
Nevertheless, least squares fits were made where we grouped Wasp Prime-9 with 
Yuma-4 and Encore-9 with Wasp-1. 

(U) There were eight non-ideal air bursts within the scaled HOB from ‘ 
about 40 to 73 metres. Various combinations of these events were used for 
obtaining least squares fits. We combined light dust events with heavy dust 

events in an attempt to evaluate any significant differences between these two 
categories. After establishing which of the events were light dust or heavy 

dust blast waves we used the Tq values derived from the asphalt surfaces for 
comparisons (see Table 4.4) with the light dust events. 

(U) A total of 23 least squares fits were made to the scaled Iq versus 

scaled ground range data. Table 4.5 lists the combination of events, the 

burst conditions, the number of data points used in the fit, the derived 

functional least squares equation, the associated statistical data, and 

comments regarding what was used in the fits and the scaled range applicable | 
to the fit. Under the remarks column we make a distinction between Tg and 
Iq measured. The Ig is, of course, the average of lay and Igavg and I, 
measured are the measured values only. The I, measured were used as addi- 
tional points for the fits as well as to establish the validity of the 
measurements. In one case we fitted the data using only Iq measured (Met-12 
and Apple II-13). The values for these events are given in Table 4.6. Also, 
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in a couple of cases we fitted the curves by eye because of lack of data 
points, e.g., Bee-6 and Wasp Prime-9. The total number of Figures prepared 


for the various combinations of events was 18 and the results are shown in Figures 


4.1 through 4.18. These curves are shown in sequence according to the average 
burst height. 


(U) Two least squares fits were made for the surface bursts (Figure 4.1). 


The reason for this is that the functional form selected did not fit the data 
points at the larger ranges. Fit 1 and Fit 2 are of the same functional form 
but Fit 1 contains all data points used and Fit 2 contains 10 of these points 
which were arbitrarily selected to fit the data points at the larger ranges. 
We place high confidence on these fits because of the high coefficient of 
determination, Re, (see Table 4.5) and the low values for the ¢ itistical 
parameters and in particular the percentage errors, E;(%), which are about 

13 percent. Again, these scaled Iq values represent the ideal/near-ideal 
blast waves. 

(U) The next set of curves, Figures 4.2 through 4.6 rep: ~ent the 
scaled lg for non-ideal blast waves. The cverage range for the HOB is from 
41.5 metres out to 44.6 metres — an arbitrary selection for the HOB. Note 
on Figure 4.2 that the asphalt data on Met-12 fits reasonably with the data 
of Turk-4, which we later defined to be light dust. The asphalt data includes 
the lgay calculated from vehicle displacements over asphalt surface using 
the quadratic fits to the Iq versus displacement which were obtained for the 
desert surface (Nevada Test Site) and for the coral surface (Pacific), 
Section (3.1.2.4). The Iq versus displacement for the asphalt surface are 
comparable to the desert and coral surfaces (see Section 3.5.3). 

(U) Figures 4.3 and 4.4 combine what we consider to be light dust and 
heavy dust-laden blast waves. Met-12 Iq values are generally higher than the 
Turk-4 or Apple II-13 Iq values (Figure 4.3). The Turk-4 data and Apple II- 
13 data appear to coalesce. Matter of fact, when we fit just the Iq values 
on Turk-4 and Apple II-13 we get good agreement (Figure 4.4, Fit 2). This 

is not the case when we included the Ig measured values (Figure 4.4, Fit 1). 
The trend here shows an increase in Ig values at closer in ranges for Apple 
II-13 than for Turk-4. We did note that as reported in Reference 9 (WT-1123) 
the damage to vehicles (likewise displacements) on Apple II-13 was not as 
extensive as for other events. The reason given was that the vehicles were 
located on a surface which was unusually soft and sandy. We proceeded to 
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compare only Apple II-13 and Met-12 wherein the surface was also unusually 
soft and sandy, particularly along the desert blast line. These comparisons 
are shown in Figures 4.5, Fits ] and 2 and Figure 4.6, Fit 3. The scatter of 
data is about the same for Fits 1 and 3 and somewhat larger for Fit 2 (see 
Table 4.2). When deriving the Iq values for Apple II-13, the Tdavg derived 
from vehicle displacements and averaged with Iq (measured) tended to lower 
the Iq values relative to Igo: The opposite of this tendency was the case 
for the majority of Iq on other events, i.e., the Ig, avg derived from vehicle 
displacements averaged with Iqo values increased Ig relative to Tag: We, 
therefore, compared only the measured Ig on Apple II-13 with Met-12, Figure 
4.6. Above Ig = 10 kPa-sec the comparison is poor but for the remaining Iq 
values there is good agreement. Note that the percentage difference in range 
for the three Fits is small for Tq between 0.8 kPa-sec and 10 kPa-sec. If we 
believe the measurements (our inclination is that they are valid and reason- 
able) then we have a dichotomy with regard to Apple II-13. 

(U) On the one hand, just using Iq values for Apple IT-13 and combining 
with Turk-4 indicates that Apple II-13 could be considered a light dust event 
or we could argue that Turk-4 is a heavy dust event. However, Turk-4 Iq data 
at close in ranges are definitely below the Met-12 data (one of the measures 
for heavy dust and light dust) so Turk-4 was classified as light dust. Apple 
II-13 when combined with the Met-12 data could be classified as a heavy dust 
event (ali previous evaluations stress the fact that Met-12 was a heavy dust 
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event) since there is some agreement when the data from these two events are 
combined. We dwell on this matter to point out that the data per se for a 
number of events does not permit a clear distinction between light dust and 
heavy dust blast waves throughout the Iq versus range curves. All one can 
do is surmise that the type of surface conditions will produce varying 
amounts of dust content in the blast wave. 

(U) Now we come to the I, versus range curves, Figures 4.7 to 4.13, 
which we classified as light dust events since all the Ig data on these 
curves are lower, particularly at the close in range than the heavy dust 
events. These curves are for HOB ranging from 46.8 metres out to 66.9 
metres. If one were to superimpose these curves on top of one another they 
would not vary by more than 13 percent in range and this is mainly in the 
region of iy from 0.6 kPa-sec to 0.8 kPa-sec. As before, in this set of 
curves we made use of the Iq measured as well as Iq. It is to be noted that 
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the asphalt results, Met-12 have also been included with Turk-4 and Smoky-15. 
Here, we classify Smoky-15 as a light dust event. The asphalt results fit 
reasonably well with the other cata (Figure 4.7). We observe the same situa- 
tion when we included asphalt results from Bee-6 with Moth-2 and Grable-10 
(Figure 4.13). Again the asphalt Tq values from Bee-6 include the ldavg from 
vehicle displacements over the asphalt using the quadratic fits (Section 
3.1.2.4). An interesting observation about Grable-10 is that this event was 
detonated in the Frenchman's Flat area as was Met-12, yet we consider Grable- 
10 as producing light dust-laden blast waves. We pointed out in our initial 
evaluation (Reference 1) that on Met-12 the vehicle response along the desert 
line (disturbed, very fine loose soil) was greater than the vehicle response 
along the Desert Rock sector (same soil but undisturbed). For Met-12, the 
Igavg (calculated from vehicle displacements) was greater on the desert line 
than on the Desert Rock sector. We can safely assume, therefore, the higher 
Idavg along the desert line was the result of the disturbed surface (producing 
higher dust concentrations in the blast wave). The Frenchman's Flat soil con- 
sists of alluvial deposits which are fine microsize particles similar to talcum 
powder. The desert line region was disturbed by construction of test items, 
digging up cable ditches, and high concentration of daily traffic prior to 
detonating the device. These are somewhat unique conditions. ‘No such acti- 
vity took place in the Desert Rock sector. The vehicles on Grable-10 were 
located in the Desert Rock sector and thus, had lower values of Iq equivalent 
to the Iq values of Moth-2, Apple I-8 and Smoky-15, which were classified as 
light dust events (see Figure 4.11). We also note that Fits 1 and 2, Figure 
4.11 nearly coincide over the total range of Tq values and this is the same 
situation for Fits 1 and 2, Figure 4.10. 

(U) The next curve, Figure 4.14, is a fit to the data obtained on 
Grable-10, Moth-2 and Bee-6 (Fit 2). Two of the if data points do not signi- 
ficantly alter the curve, however, as for example, Fit 1 for Moth-2 and Bee-6, 
Figure 4.15. The curves of Fit 2 (Figure 4.14) and Fit 1 (Figure 4.15) are 
about the same. Since these two points did not necessarily influence the fit 
we decided to fit the data points by eye. The reason for this is that we 
assumed that Bee-6 could be considered to produce heavy dust-laden blast 
waves. The desert surface on Bee-6 was a fine, soft sandy soil similar to 
Apple II-13. However, the four data points on Bee-6 do not necessarily 
validate this assumption. 
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(U) The next set of curves, Figures 4.16, 4.17 and 4.18 deal with near- 
ideal blast waves. The average HOB for these curves range from 127 metres up 
to 240 inetres. The HOB = 127 metres is represented by Yuma-4 (HOB = 109 
metres), and Wasp Prime-9 (HOB = 145 metres) shown in Figure 4.16, Fit 1. We 
do not feel that this curve, Fit 1, is representative for the Iq values at 
this average HOB for several reasons. First, we mentioned previously that we 
questioned the I, values on Yuma-4; second, the scatter of the data is large 
between those two events; and third, the large difference in HOB (36 metres) 
between these two events could be a factor in the noted scatter. We feel, 
therefore, that in spite of the lack of data, the fit by eye to the Wasp 
Prime-9 event is more reasonable. The remaining two curves, the combined 
events, Encore-9 and Wasp-1, Figure 4.17 and Encore-9, Figure 4.18, do not 
show any substantial difference. The maximum difference between these curves 
is less than 9 percent and that falls close to the scatter for each curve 
(see Table 4.2). Because of the uncertainty of the data on Wasp-1 (vehicle 
orientations and smal] displacements which poorly determined Igayg values) we 
consider the Erncore-9 curve only to be more reliable. 

(U) To recap, we place high confidence on the Iq data for the surface 
bursts (ideal/near-ideal) and for the HOB range from about 40 to 70 metres 
(light dust). There is less confidence in the height of burst at HOB from 
109 metres out to 240 metres. But this is only because of the small amount 
of data available at these burst heights. Nevertheless, the Iq values seem 
reasonable, except for Yuma-4. The heavy dust events, represented by Bee-6, 
Met-12 and Apple II-13, indicate Iq values to be 2 to 3 times the values for 
the light dust events (compare Figure 4.5, Fit 1, Heavy Dust with Figures 
4.2 and 4.7). The surface conditions for the heavy dust events producing 
these blast environments were unique, particularly the Frenchman's Flat area. 
Only under these conditions would a similar blast environment be produced and 
these surface conditions would be rare to find in many other parts of the 
world. We tend to believe the validity of the measurements (there are some 

exceptions of course) but the Iq which combine measurements and laavg from 
vehicle displacements should, in general, be more reliable which is what we 


wanted to show. 


4.2 (U) DYNAMIC PRESSURE IMPULSE, Iq: HOB CHARTS 
(U) We presented various scaled Iq versus scaled range curves for three 
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types of blast environments; ideal/near-ideal, light dust, and heavy dust 
(see Section 4.1). The procedures for classifying each event or a group of 
events into one of the above categories were outlined. Use was made of these 
curves for each type of blast category to prepare Tables which list the I 
values as a function of range for different burst heights. The scaled I 
values in these Tables extend from 0.4 kPa-sec out to 10 kPa-sec. The scaled 
range for each Ig value was selected by eye from the fitted curves. Tables 
4.7, 4.8, and 4.9 give this information for ideal/near-ideal, light dust and 
heavy dust blast environments, respectively. The Tp values given are for a 
height of 0.91 metres (3 feet) above the ground surface. 

(U) The Iq versus range data for the ideal/near-ideal blast waves are 
given for five burst heights (Table 4.7). Two of these burst heights include 
the Iq values from Yuma-4 (HOB = 127 metres) and Wasp-1 (HOB = 224 metres) 
which yielded questionable data (see Section 4.1). However, a comparison of 
the ranges for several Ig values for HOB = 127 metres and 145 metres shows 

the percentage difference is less than 9 percent. The exceptions are for Tq 
of 5 kPa-sec and 4 kPa-sec. Here, the differences are about 21 percent for 

Tq = 5 kPa-sec and 14 percent for Ig = 4 kPa-sec.. 

(U) The next set of Iq versus range data is for the light dust events 
(Table 4.8). The averaged HOB in this set of data extends from 41.5 metres 
out to 66.9 metres. If we neglect the extrapolated values of Iq = 10 kPa-sec, 
the maximum percentage difference in range for given values of Iq is about 17 
percent for these burst heights and this difference is due primarily to the 
extrapolated values for the lower if values. If we further neglected all of 
the extrapolated values then the percentage difference would be about 10 
percent. This indicates that these data could be consolidated into one aver- 
aged burst height and least squares fitted without introducing large errors. 
(U) The heavy dust Tq versus range data spans the averaged HOB from | 
43.4 out to 73.2 metres (Table 4.9). In this set of data, there are some 
events that we classified as light dust which have been combined with the 
heavy dust events. In Section 4.1, we indicated the difficulties encountered 
in attempting to distinguish between light dust and heavy dust events over 
the total Ig versus range curves. The events which we feel clearly fall into 
the heavy dust category based on the measures used for distinguishing light 
and heavy dust, are Met-12 and Apple II-13. The surfaces for both these 
events were fine, loose sandy soil and the magnitudes of dynamic pressure 
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impulse were 2 to 3 times greater than for the light dust events. Bee-6 
could fall into this category as well. Similar surface conditions existed 
but the maximum difference in Iq was only about 1-1/2 times greater than for 
light dust and this was over a small portion of the range (compare Bee-6 and 
Moth-2, Figure 4.15). The set of Tq versus range data from Table 4.9 that 
could then be interpreted to be representative of heavy dust events are only 
two, HOS = 44.6 metres and 73.2 metres. Between these two burst heights the 
difference in range for constant Iq is small, less than 8 percent. Even over 
this whole span of HOBs, neglecting the extrapolated values, the range differ- 
ences for constant Iq values decreases from about 29 percent down to 10 
percent with decreasing Iq. 

(U) For constructing the HOB Charts, the Iq versus range data from the — 
three types of blast conditions have been combined and placed into two cate- 
gories — light dust/near-ideal and heavy dust/near-ideal blast waves. The 
HOB Charts are shown in Figures 4.19 and 4.20, one for light dust/near-ideal 


{v. and the other for heavy dust/near-ideal blast environments, respectively. 
‘2 The burst heights selected for constructing these charts were: three for the 
oy near-ideal, 0, 145, and 240 metres; two for the light dust, 49.2 and 61.1 


metres (Fit 1); and two for the heavy dust, 44.6 (Fit 1) and 73.2 metres (eye- 
fit). The selection of only two HOBs for the light dust events was mainly to 
minimize clutter of data on the charts. As we pointed out for the light dust 
events there were small differences in Iq versus range curves over the span 

of HOB from 41.5 to 66.9 metres. As mentioned previously, only two of the 

Iq versus range curves could be interpreted to be heavy dust events. Thus, 
the five burst heights selected from each combination were sufficient to, at 
least, determine the trend of the iso-Ig contours. The error bands have not 


been placed in the contours but these can be readily obtained from the statis- 
tical data derived by least squares fit procedures (see Table 4.5). 

(U) Upon inspection of these Charts we note the effect of HOB on range 
for equivalent i values. As the burst height increases, the range initially 
increases and then decreases, but for some tow values of Iq the range 
increases with increasing HOB. Moreover, the iso-I, are not everywhere para- 
llel. We also note that from about HOB = 40 metres to 70 metres the range is 
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(U) Although the Iq values from the heavy dust events (HOB = 44.6 

metres) were 2 to 3 times higher than the light dust events (HOB = 49.2 

metres) the comparison of Tq as a function of range shows that the range 

difference is about 27 percent greater for the heavy dust events at the 

higher Tq than for the light dust events. The range difference then 

decreases to about 10 percent at the lower Iq. As the burst height is 

increased up to about 70 metres the range difference between these two types 

of blast environments becomes smaller, down to about 6 percent. These curves, 

heavy dust and light dust, converge at the lower Tq values. 

(U) Here, we are comparing the heavy dust data from three events which 

are unique and one of these (Bee-6) has only limited data with data from five 
light dust events which we feel are reasonably reliable. Under these circum- 
stances and since the range differences are not that large, our tendency is to 

place more importance on the HOB Chart which includes the light dust events 

than the chart which includes the heavy dust events. Other observations, but 
somewhat qualitative, which further lead us toward this trend are the surface 

conditions of the light dust events, the pressure-time measurements, and the 
agreement between Ty derived from the light dust events and the asphalt sur- 
faces of Bee-vu wid Met-12. The light dust events were detonated over surfaces 

which consic ed of rough gravel mixed in with sand and, in some cases, vegeta- 
tion was interspersed throughout the area. For all of these events except 
Smoky-15, only BRL mechanical self-recording gages were used along the blast 
Tine. Cable ditches or any major obstruction were not required to install 

the gages. The surface conditions were more or less undisturbed. Additionally, 
a cursory examination of the pressure-time measurements made on the light dust 
events shows less high frequency oscillations than for the heavy dust everts. 
This implies that the influence of dust on these measurements (light dust) 

was minimal. We, therefore, believe that the main influences on the derived 

Ig for the light dust events was the result of precursor action and not | 
necessarily the dust. This is somewhat borne out by the [4 values obtained 

over the asphalt surfaces on Bee-6 and Met-12. The asphalt Tq were combined 

with Ig from the light dust events and both are represented reasonably well 

by a single curve for that specific averaged burst height (see Figures 4.2 

and 4.13). No significant amount of dust was present over the asphalt sur- 

face. For these reasons we tend to consider the HOB Chart which includes the 
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light dust to be more representative of the expected blast conditions over 
various surfaces where precursors (non-ideal) would be formed. 

(U) Since the surface conditions (disturbed, micro-size particles) of 
the heavy dust events such as on Frenchman's Flat area were unique, the blast 
environments produced were also unique. These unique surface conditions gave 
rise to heavy dust loading of the blast wave (precursor) producing much higher 
Tg than for the light dust events. We, therefore, tend to place less impor- 
tance on the heavy dust events and perhaps even to disregard them. 
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TABLE 4.1 (U) Dynamic Pressure Impulse versus Range, 
Final Averages -- Ideal/Near-Ideal. Blast Waves, Surface Bursts 


UNCLASSIFIED 
ACTUAL SCALED I ACTUAL SCALED 
RANGE RANGE lag avg Iq Iq 
OPERATION/EVENT Metres Metres kPa-sec kPa-sec kPa-sec kPa-sec 
Castle, Koon-3 198] 372.2 2.83 2.29 2.56 0.494 
2499 469.6 1.26 1.62 1.43 0.276 
3322 624.2 0.647 -- 0.647 0.125 
4206 790.3 0.346 -- 0.346 0.067 
Castle, Nectar-6 2316 193.6 33.9 22.1 28.0 2.4) 
2408 271.3 29.4 21.5 25.4 2.18 
2591 216.6 24.7 21.8 23.0 1.96 
2743 229.3 20.7 19.8 20.3 1.74 
2920 244.1 17.7 -* 7.7 1.92 , 
3990 333.6 7.73 §.38 | 6.56 0.563 
Redwing, 762 223.3 6.18 8.05 7.12 2.14 
Lacrosse-] 844 247.3 4.78 5.58 5.18 1.56 
1021 299.2 2.89 2.92 2.90 0.873 
1189 348.4 1.82 2.22 2.02 0.608 
1334 390.9 1.34 1.49 1.42 0.427 
Redwing, Zuni-3 2135 140.3 77.4 62.5 70.0 4.72 
2530 166.2 51.4 30.4 40.9 2.76 
3170 208.3 28.3 37.4 32.8 2.21 
3566 234.3 19.8 29.5 24.6 1.66 
4206 276.3 1253 14.6 13.4 0.903 
5029 330.4 6.71 4.79 §.75 0.389 
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TABLE 4.2 (U) Dynamic Pressure Impulse versus Range, 
Final Averages -- Ideal/Near-Ideal Blast Waves/Air Bursts 


CONFIDENTIAL - FORMERLY RESTRICTED DATA 


ACTUAL SCALED ACTUAL SCALED 


RANGE RANGE a9 Mavg Tq Tq 

OPERATION/EVENT Metres Meters kPa-sec  kPa-sec kPa-sec kPa-sec 
Upshot-Knothole, 756 245.2 3.52 3.17 3.38 1.24 
Encore-9 872 282.9 3.00 2.84 2.92 1.07 
933 302.7 2.75 3.01 2.88 1.06 

“198 388.6 1.82 2.45 2.13 0.781 
1329 431.1 1.51 1.64 1.58 0.579 

1704 552.8 0.830 0.870 0.850 0.312 

1996 647.5 0.680 0.759 0.719 0.264 

Teapot, Wasp-1] 354 317.3 0.917 0.979 0.948 0.978 
393 352.2 0,779 Meld 0.951 Q.981 

434.5 389.4 0.643 0.80 0.722 0.745 

460.8 413.0 0.567 1.05 0.808 0.834 
607.3 544.3 0.315 0.688 0.502 0.518 

Teapot, Wasp 237 152.3 8.63 9.24 &.94 6.83 
Prime-9 311 199.8 6.75 8.38 7.57 5.78 
465 298.8 1.82 -- 1.82 1.39 

618 397.1 0.871 -- 0.871 0.665 

Redwing, Yuma-4* 76 132.8 2.38 3.23 2.81 5.07 
107 187.0 1,30 1.85 1.56 2.82 

122 213.3 2.60 


1.63 1.85 1.44 


(Uy *Yield from the above information makes this event CONFIDENTIAL - 
FORMERLY RESTRICTED DATA 
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TABLE 4.3 (U) Dynamic Pressure Impulse versus Range, 
Final Averages -- Non-Ideal Blast Waves/Air Bursts 


UNCLASSIFIED 
ACTUAL SCALED I , ACTUAL SCALED 
OPERATIGN/EVENT Metres Metres kPa~sec kPa-sec kPa-sec kPa-sec 
Upshot-Knothole, 736 287.7 5.13 5.66 5.39 2.37 
Grable-10 844 329.9 2.69 2.53 2.61 1.11 
(Desert Rock) 1335 433 ae 0.686 0.686 0.302 
Teapot, Moth-2 287 204.1 6.75 = 6.75 5.39 
366 260.3 3.35 3.86 3.6] 2.88 
4}} 292.3 2.42 2.05 2.24 1.79 
457 325.0 1.88 1.47 1.68 1.34 
518 368.4 1.35 a 1.35 1.08 
Teapot, Turk-4 716 193.0 20.7 am 20.7 6.51 
914 246.4 11.8 10.5 11.2 3.52 
1030 277.7 7.69 9.81 8.75 2.75 
1128 304.1 5.66 7.77 6.72 2.11 
1280 345.1 2.32 «- 2.32 0.730 
Teapot, Bee-6 549 263.7 8.10 9.36 8.73 4.79 
(Desert) 610 293.0 4.17 4.44 4.3) 2.37 
701 336.7 2.32 2.59 2.46 1.35 
777 373.2 1.91 2.66 2.28 1.25 
Teapot, Aople I-8 527 205.6 17.3 i. 17.3. 7.95 
622 242.6 9.10 -- 9.10 4.18 
902 351.9 1.95 2.62 2.29 1.05 
99] 386.6 1.20 1.77 1.48 0.68 
1128 440 0.38 0.76 0.57 0.262 
Teapot, Met-12 610 208.9 43.0 41.0 42.0 16.4 
(Desert) 686 235.0 37.6 41.1 39.3 15.3 
762 261.0 16.8 15.7 16.3 6.36 
838 287.0 13.8 15.2 14.5 5.66 
914 313.0 3.51 3.57 3.54 1.38 
Teapot, Apple II-13 518 160.3 41.9 = 41.9 15.4 
625 193.4 34.0 -- 34.0 12.5 
a) 808 250.0 19.3 “= 19.3 7.08 
hey 914 282.8 10.3 7.81 9.06 3.32 
ave 1006 311.3 5.54 3.66 4.60 1.19 
Z ‘ 1128 349.0 3.98 3.23 3.61 1.32 
of 1219 277.2 2.55 2.81 2.68 0.983 
on Plumbbob, Smoky-15 457 122.6 138 Be 138 44.7 
726 194.8 28.5 36.8 32.6 10.6 
* 789 211.7 23.7 -_ 23.7 7.68 
ne 84] 225.6 16.0 14.1 15.1 4.90 
at 897 240.7 12.5 7.70 10.1 $,27 
1038 278.5 6.90 7.22 7.06 2.29 
1181 316.9 3.75 3.59 3.67 1.19 
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TABLE 4.4 (U) Dynamic Pressure Impulse versus Range, 
Final Averages, Asphalt Surface -- Non-Ideal Blast Waves/Air Burst 


UNCLASSIFIED 


ACTUAL SCALED I I ACTUAL SCALED 
RANGE — RANGE Tt) Favg 


I 
OPERATION/EVENT Metres Metres kPa-sec kPa-sec ipalsec kPadsec 


al a a el 


Teapot, Bee-6 549 263.7 6.59 6.97 6.78 3.72 
610 293.0 4.16 2.51 3.38 1.85 
701 336.7 1.60 2.10 1.85 1.02 
777 373.2 1.02 1.24 1.13 0.62 | 
Teapot, Met-12 — 610 208.9 19.0 18.9 19.0 7.42 | 
686 235.0 13.7 14.8 14.2 5.54 
762 261.0 8.40 8.34 8.37 3.27 | 
338 287.0 5.58 4.38 4.98 1.94 
914 313.0 Vive 1.37 1.54 0.60 | 
| 
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TABLE 4.6 (U) Dynamic Pressure Impulse versus Range 
for Met-12 and Appie II-13 -- Measured Values 


UNCLASSIFIED 
ACTUAL I SCALED I 
RANGE qo RANGE % 
Metres kPa-sec Metres kPa-sec 
MET-12 

610 48.0 208.9 18.7 
686 37.2 235.0 14.5 
762 17.2 261.0 6.71 
838 17.5 287.0 4.49 
914 3.4 313.0 1.33 

APPLE II-13 
518 41.9 160.3 15.4 
625 34.0 193.4 12.5 
808 19.3 250.0 7.08 
914 12.1 282.8 4.44 
1006 7.4 311.3 2.71 
1128 4.8 349.0 1.76 
1219 2.48 377.2 0.910 
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TABLE 4.7 (U) Height-Of-Burst For Dynamic Pressure Impulse Versus 
Ground Range At 0.91 Metres Above Surface 


Ideal /Near-Ideal 
UNCLASSIFIED 


verage HOB tres 
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* Extrapolated 


123 


UNCLASSIFIED 


ei HE 4 Se Ba 2 ee ee ea eee 


eee er 


er 
- 
“5 
wy 
be 
pice 
- ee are 
rae) 
poze Lodeaqxj, Shs 
SLY 6LY BLY - diy 9lp +l0b s€8€ «BSE v0 Aes 
3 06€ ese ose BLE LE | ESE EME 9°0 A 
S9E oe OE BSE Sse e0SE—«OEESCi«CEE 8°0 £% 
rE ySe «I: dla 8D lve [ld 8D. Lge BEE = «SEE 2ESC—*«E o'l us, 
a vee Ove ous Sze aug B2E Eze 40LzE—iéiE:C:*«SSCLE 27k a 705 
w lie ee Lte “AEC (it OSC(tiOESC(‘ié‘iZ:*SC<“‘«é‘ét SCS] = as 
x £62 662 qnoqy 062 qnoqy 162 062 882 282 882 2 2 oe 
a 092 i9z 292 992 e972 09-——isiBSZ—(tsCSZ € a “25 
a Sez nd Sez “Geo Sb2ZCLZSséZ b 3) Bs 
> slz 012 lez lez fz ZZ RACV S rs a8 
= «ISL -- 061 261 O6L ssl O6L +591 ol *< 
ee 
2 itd @ 41d i tJ 2 Lid tis rE) ne 
6°99 L*€9 ig U'tg +85 7 a a he) Ae I eX 
$8430) GOH abeusAy 5 
Q314ISS¥TONN ae 


{eapl-uon/ysng 74614 


BIVpANS BAOGY SA1J9W 16°0 Ty ebuey puNnouyg snssaq 
as_ndwy aanssaag Dpweudg 204 354Ng-JQ-24619H (N) BF JTL 


ioe 
Shp 


WY 


—_ 


UNCLASSIFIED 


414~-3A3 
cel 


= ver -- *00P *O0b¢ += == 
== S6E == #18€ *S0¥ ne =< 
ra Ze SLE 89€ +#l8€ ¥26E 9L€ 
L8€ LSE S9E 8SE ¥9E €Le 6SE 
o9¢e eve LSE Ose ose LS€ StE 
cee Lee Lve OvE Ste LEE 62E 
OO€ €0€ LEE GZE Lle OLE LOE 
292 892 ole £0€ 62 4 Lle 
SEZ Ove 262 882 622 Ove 9S2 
SIz SLz 8l2 2l2 692 ar 8Ez 
< #21 822 S22 8E2 =- 061 

L Lid 2 itd € 14 2 iid t 11d 2 Ltd t 114 
1°69 6°99 9° tb 9°tP 9° bP L°vb L* vt 


S9439W) GOH abeuaay 


G314ISSYIONN 
Leap]-uoN/3sng Aaeay 


BICJANS BACT’ SA1ZW 16°0 3Y abuey punouy snssaa 
asindw] ainssaug J,weukg soy 4S4ng-40-2461aH (A) 6°b 3JTGVL 


paze|odeazx3, 
es ¢°0 
== 9°0 
#L6E 8°0 
69€ o'l 
BPE Zl 
O€E S*t 
80€ 2 
$82 € 
022 v 
092 S 
LEZ Ol 
9S-egy 
yep by 


125 


UNCLASSIFIED 


: BAWMAeMRA RA BREF RRB en SEee ZARA RR BRA AS AB Rm OE em) me ee a es | oS tt Ow ee ee ewe EER RRR MMAR PR AANA RS! ei Wt 


UNCLASSIFIED. 
'Or . the Bass ert ? = « é = See ee 
po IT Tt Tg 3.24 - 0.0137x + 0.000008 x2 
AN ee Se et Re 2 0.980 
FIT 2: £n Ig = 1.945 - 0.0085x + 0.000034 x? 
R? = 0.961 
Zs “~~~ # Castle, Koon-3 
“~~ - @ Castle, Nectar-6 
3 w Redwing, Lacrosse-] 
7 ane © Redwing, Zuni-3 
2 (ace Semeiaeps Sr tema id Aan en naam Sarre ions til Slows Gee ee 9 
10° enon 
e : — ven ens Dee, 
——s —— 
“ . 
- 4 
2 
a 
‘ - 
Vee —_ 
= | 
ev 
. 2 
aa, 
1071 
3 
5 pS 
3 os ee 2 ee - ame aie rc, eee 
1072 Sealed Range - Metres 
100. 200. 4300  # «400 506° ° ~~ 600 ° © ~~ 700 300 


UNCLASSIFIED 


FIGURE 4.1 (U) Dynamic Pressure Impulse vs. Range - Surface Bursts 
Scaled to 1 KT 
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FIGURE 4.2 (U) Dynamic Pressure Impulse vs. Range -- Average HOB = 
41.5 Metres, Scaled to 1 KT 
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FIGURE 4.3 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 
43.4 Metres, Scaled to 1 KT 
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FIGURE 4.5 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 
44.6 Metres, Scaled to 1 KT 
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FIGURE 4.6 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 
44.6 Metres, Scaled to 1 KT 
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FIGURE 4.7 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 46.8 
Metres, Scaled to 1 KT 
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FIGURE 4.8 ‘U) Dynamic Pressure Impulse vs. Range - Average HOB = 49,2 
Metres, Scaled to 1 KT 
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FIGURE 4.9 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 52.6 
Metres, Scaled to 1 KT 
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FIGURE 4.10 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 58.4 
Metres, Scaled to 1 KT 
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FIGURE 4.11 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 61 
Metres, Scaled to 1 KT 
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Figure 4.12 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 63.7 
Metres, Scaled to 1 KT 
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FIGURE 4.14 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 66.9 
Metres, Scaled to 1 KT 
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FIGURE 4.15 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 69.1 
Metres, Scaled to 1] KT 
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FIGURE 4.18 (U) Dynamic Pressure Impulse vs. Range - Average HOB = 240.0 
Metres, Scaled to 1 KT 
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SECTION 5 
(This Section is Unclassified) 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


5.1 SUMMARY AND CONCLUSIONS 


5.1.1 Dynamic Pressure Impulse versus Displacement 
Displacement versus Dynamic Pressure Impulse 


Two types of functions have been used for least squares fit to 
dynamic pressure impulse versus displacement and vice versa. The power func- 
tions were used to separately fit the data for the WWII and the M38Al 1/4 ton 
trucks. After further analysis in which we pointed out that these two vehicles 
could be treated as one type of truck, we combined the data for the 1/4 ton 
trucks as well as for the 2-1/2 ton trucks (M35 and M135) and used quadratic 
functions for least squares fits. Also, the Iq values for the non-ideal 
events were iterated and used in the latter least squares fits. The resultant 
equations for Ig versus displacement are: 

Combined Vehicles: M38A1 and WWII - 
SO fn Iq = 0.133 + 0.376 £n D + 0.0566 (£n D)? 
FO ln Iq = 0.732 + 0.465 £n D + 0.0220 (£n D)® 


Combined Vehicles: REO M35 and GMC M135 - 
SO fn Iq = 0.429 + 0.333 £n D + 0.0275 (Ln D)* 


and for displacement versus Ig the equations are: 
Combined Vehicles: M38A1 and WWII - , 
SO gn D = -0.223 + 2.222 fn Iq - 0.214 (Ln Iq) 


* FO en D = -1.362 + 2.081 Ln Iq - 0.0998 (fn Iq)? 
is Combined Vehicles: REO M35 and GMC M135 - 
in SO fn D = -1.365 + 3.339 £n Iq - 0.4259 (£n Iq)* 


We find good correspondence between dynamic pressure impulse and dis- 
placement, in spite of the large scatter of data. A number of uncontrolled 
variables present in the experiments, which we described, contributed to the 
scatter. These variables are statistically averaged out in the use of least 


rote 


squares fit. One, therefore, can calculate with reasonable confidence the Ig 
from known vehicle displacements and calculate displacements (such as onset 
of motion) for given values of dynamic pressure impulse. 
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5.1.2 Displacement versus Scaled Ground Range 

For different combinations of events we fitted the displacement 
versus scaled ground range using various values of n (the exponent of W used 
in scaling ground range). The best value of n is somewhat dependent on the 
burst height and surface conditions. The values of n for best fit of data 
were from n = 0.4 to n = 0.43. Our reason for these fits was to determine 
the feasibility of scaling the ground range for displacement since we correla- 
ted vehicle damage with displacement (Reference 17). Thus, for known displace- 
ments as a function of scaled ground range damage can be then predicted for 
different yields. The utility of this scheme needs further evaluation. 

5.1.3 Dynamic Pressure Impulse (Unscaled and Scaled) versus Scaled 

Ground Range 

We again fitted the unscaled dynamic pressure impulse versus scaled 
ground range using different values of n. The values of n which best fitted 
the data were n = 0.42 or n = 0.43. These were then compared to the scaled 
dynamic pressure impulse versus scaled ground range. We concluded from an 
overall] viewpoint there is little difference in reliability between the two 
sets of curve fits to the data. The scaled or unscaled impulse versus scaled 
ground range curve fits may be used with equal confidence. 

5.1.4 Surface Effects on Dynamic Pressure Impulse, Displacement and 

Damage 

Comparing the data from water, asphalt, and Desert Rock for the 
same events, (Bee-6 and Met-12) we conclude the following: 

1) The Iq for asphalt are lower than those for water and both are 
lower than for the desert line, (Met-12). On Bee-6, again, the Ig for asphalt 
are lower than for the desert line. 

2) Displacements of vehicles on the water line are similar to 
the desert line, (Met-12) while the displacements on asphalt are smaller than 
on desert surface. The vehicle displacements on the Desert Rock are about the 
same as on the asphalt surface (Bee-6 and Met-12). 

3) Comparison of displacement versus dynamic pressure impulse 
for desert and non-desert cases shows that there is no substantial difference 
among desert, water and asphalt surfaces. 

4) Damage on the water surface was somewhat lower than on the 
desert surface while the damage on the asphalt surface on the average 
was somewhat higher than on the desert surface. The damage on Desert Rock 
147 
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for a given displacement was higher than that on the desert line. Thus, 

fer a given Iq value, displacement is at least approximately independent of 
the surface on which the vehicle is placed but for a given displacement damage 
is greater on the harder surfaces. 


We used the once-iterated Iq values (ideal/near-ideal, blast 
waves) and twice-iterated Ig values (non-ideal blast waves) for constructing 
scaled Ig versus range for several combinations of yield and burst heights. 
For most of the events we were able to distinguish the ideal/near-ideal blast 
environments from non-ideal blast environments. The non-ideal blast environ- 
ments were further divided into light dust and heavy dust events. 

We have high confidence in the scaled Iq versus range data for the 
surface bursts (ideal/near-ideal) and for the HOB span from about 40 to 70 
metres (light dust events). There is less confidence on the Ig data above 70 
metres but this is only because of the smal] amount of data available above 
this burst height. 

We concluded that the surface conditions for the heavy dust events 
were unique and similar surface conditions would be rare in other parts of the. 
world., We also concluded that dynamic pressure measurements were reasonable. 
More importantly, combining the Iq results from gages and those calculated 
from vehicle displacement results in greater reliability and accuracy in the 
Iq values. We also conclude that the Iq data obtained on the asphalt surface 
(Bee-6 and Met-12) fit reasonably wel] with data from the light dust events 
having similar burst heights. 


5.1.5 Scaled Dynamic Pressure versus Scaled Ground Range 


Two HOB charts have been constructed, one for light dust/near- 
ideal and one for heavy dust/near-ideal. The light dust and heavy dust portions 
of the charts apply to HOB from 40 to 70 metres and the near-ideal apply to 
surface bursts and to HOB above 70 metres. 

A comparison of heavy dust and light dust events shows that the 
Iq values for heavy dust are from 2 to 3 times the Iq values for light dust 
events. In terms of range, however, for equal Iq values the heavy dust is 
greater by about 27 percent at the higher Iq than for light dust and then 
they both converge to about the same range for lower Iq values. 

In the HOB region from 40 to 70 metres the range for iso-Ig 
seem to be about the same for either the light dust or heavy dust events. 
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Secause of unique surface conditions for the heavy dust events 
the blast environments produced were unique. Our tendency is to place less 
importance on the heavy dust events and perhaps even to disregard them. We 
consider the HOB charts which include the light dust events to be more repre- 
sentative of the expected blast conditions over various surfaces where pre- 
cursors (non-ideal) would be formed. We feel that the resultant Iq values for 
these events (light dust) were caused principally by the precursor action 
(based on comparison of asphalt and light dust Ig values). The conclusion 
is that the light dust/near-ideal HOB chart (nomenclature here may now be a 
misnomer) is the preferred chart for application to other targets. 


5.2 RECOMMENDATIONS 

1) A procedure similar to that presented should be used for other drag- 
type targets to extend the region for Ig to higher values than 10 kPa-sec 
and for Iq to lower values than 0.4 kPa-sec. To obtain the higher values, 
heavier and harder drag-type targets such as tanks & howitzers can be used. 
For the lower values, lighter and softer targets such as antennas & communica- 
tion equipment can be used. Additionally, different drag-type targets were 
exposed to nuclear blast at burst heights other than those used here. 

The evaluation of this information would allow greater confidence in the 
construction of the iso-Ig contours. 

2) There is a need to develop predictions of dynamic pressure wave- 
forms for the non-ideal blast environment. Current waveform predictions are 
based on the data obtained on heavy dust events which we consider to be 
unique. Application of waveform predictions is needed not only for deter- 
mining loading and response of drag-type targets but diffraction type targets 
as well. 

3) After a thoregh evaluation of data for increasing the confidence in 
constructing HOB charts for Tq there is a need to construct composite HOB 


charts representing overpressure, dynamic pressure, static pressure impulse 
and dynamic pressure impulse. Various types of targets are sensitive to 
different blast mechanisms causing failure. For targeting and vulnerability 
evaluations an HOB chart of this nature would show the advantages or disad- 
vantages against an array of targets, particularly when constructed by colla- 
teral damage (non-damage or limited damage to friendly targets). 
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APPENDIX A 
(This Section is Unclassified) 


IDEAL/NEAR-IDEAL DYNAMIC PRESSURE IMPULSE 


The procedure selected to provide us with initial dynamic pressure impulse 
values for ideal/near-ideal blast environments to be correlated with vehicle 
displacements was based on theoretical calculations and analytical techniques 
along with the use of measured blast parameters (Reference 1). From the theo- 
retical calculations and the analytical techniques we obtained the ratio of 
static pressure impulse to dynamic pressure impulse, Ip/Iqs as a function of 
overpressure, P,. The theoretical calculations consisted of the HULL code 
computations for Dice Throw ard for 1 KT nuclear detonation at HOB = 60 yi/3 
The analytical technique used was that of Ethridge (Reference 4) for a 500 
ton HE hemispherical charge. These ratios of Ip/Iq as a function of over- 
pressure are shown in Figure A.1. 

We pointed out in Reference 1 that the best measurements of blast para- 
meters which were reasonably reliable were overpressure and static pressure 
impulse. Thus, for each ideal/near-ideal blast event, we obtained the average 

“values of the measured overpressure and static pressure impulse at each vehi- 
cle location. These values were then used in conjunction with Figure A.1 to 
obtain the dynamic pressure impulse. 

The curve designated as 1 KT 60 yw/s (Figure A.1) was used to derive the 
ratios of Ip/Tq as a function of P, for Upshot/Knothole, Encore-9 and Redwing, | 
Yuma-4. The scaled height of burst for these two events were: Encore-9, 

240 W V/s metres and Yuma-4, 109 yi/3 metres. The curve, NHE, was used for 
Castle, Koon-3 and Nectar-6; Redwing, Lacrosse-1] and Zuni-3; and for the 

Canadian 20 Ton and 100 Ton events. Both of the Canadian events were hemi- 

spherical and burst at ground surface. The nuclear évents were surface bursts 

and it is to be noted that the NHE curve is a good approximation of a 1 KT 

nuclear weapon burst on the surface. And, of course, for Dice Throw the HULL 


calculations were used. 

The values of dynamic pressure impulse derived by this procedure were 
checked using the Rankine-Hugoniot and empirical relations. The percentage 
difference between these two procedures in the majority of cases was about 10 
percent, although some were as large as 25 percent. 

The initial estimates of dynamic pressure impulse for the nine events 
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are shown in Figures A.2 through A.10 (solid lines). We also plotted the 
first iteration of Tq on these Figures (see Section 4.1). We note again that 
Taayg which was combined with lag to obtain Iq, was derived from vehicle dis- 
placements using the quadratic fits (Section 3.1.2.4). Overall the agreement 
between the initial estimates and Iq are good with some exceptions. One of 
the exceptions is Yuma-4 which we discussed previously (Section 4.1). Other 
exceptions are 1 or 2 points such as on Nectar-6 and Zuni-3 where the Iq is 
+ 18 to 22 percent from the initial estimates. 
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FIGURE A.1 (U) Ratio of Ip/Tq versus Overpressure 
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FIGURE A.2 (U) Dynamic Pressure Impulse versus Range - Upshot/Knothole, 
Encore-9 
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—— Initial Estimates 
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FIGURE A.3 (U) Dynamic Pressure Impuise versus Range - Castle, Koon-3 
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we FIGURE A.4 (U) Dynamic Pressure Impulse versus Range - Castle, Nectar-6 
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REDWING, Lacrosse-1 
—— Initial Estimates 


r) ig - Averaga (Displacement & 
Initial Estimates 


RANGE theties) 


800 1000 1200 1400 : 
UNCLASSIFIED 
FIGURE A.5 (U) Dyiamic Pressure Impulse versus Range - Redwing, 
Lacrosse-1 
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FIGURE A.6 (U) Dynamic Pressure Impulse versus Range - Redwing, Zuni-3 
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FIGURE A.7 (U) Dynamic Pressure Impulse versus Range - Redwing, Yuma-4 
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FIGURE A.8 (U) Dynamic Pressure Impulse versus Range - Canadian 20 Jon 
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FIGURE A.9 (U) Dynamic Pressure Impulse versus Range - Canadian 100 Ton 
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FIGURE A.10 (U) Dynamic Pressure Impulse versus Range - Dice Throw 
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APPENDIX B 
(This Section is Unclassified) 


DETERMINATION OF DYNAMIC PRESSURE IMPULSE FOR NON-IDEAL SHOTS 


B.1 DATA REDUCTION PROCESS m 
The shots involved in this section were all part of Operation Teapot. 
We used two basic procedures and some variants of these to determine dynamic 
pressure impulse from the measured pressure-time histories. The measurements 
consisted of total head overpressure and static overpressure versus time, the 
latter at two locations: one at ground level and one at the same height, 
usually 0.91 metres (3 feet), as used to measure total head pressure. In 
addition, previously conducted reduction of this data (using only the two 
gages at the same height, not the ground level measurement) provided dynamic 
pressure versus time waveforms (not corrected to account for the perturbation to 
the pressure undergoing measurement by the gage used to make the measurement). 
In the following we describe the correction process for the dynamic pressure 
provided and an independent data reduction method we used to calculate dynamic 
pressure impulse from the original gage data. 
We use the following nomenci ature (Reference 3): 
Pe - absolute free stream static pressure 
Pr - absolute free stream total pressure 
Py - absolute total head Pitot pressure 
Py - ambient pressure 
APs = Ps - Pg — free stream static overpressure 
APT = Pr - Po — free stream total overpressure 
APy = Py - Po — total head overpressure 
q - dynamic pressure 


Q; - dynamic pressure impulse; subscript i denotes ith method of 
obtaining result (4 values of i used) | 


y - ratio of specific heat at constant pressure to specific heat 
at constant volume for air 


M - local free stream Mach No. of flow behind blast front 
At; - ith time interval (in a pressure-time waveform) 
G(M) - Mach correction factor 

G(M) - average value of G(M) over pressure-time history 


® There was one exception to this: Plumbbob, Smoky-15. 
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W; - weighting factor for ith time interval 


To determine dynamic pressure impulse when q(t) is provided, without 
Mach correction, we must calculate the Mach correction and integrate the 
corrected dynamic pressure over time. This involves some of the operations 
used in calculating the dynamic pressure impulse from the raw data (pressure 
waveforms), so we start with the latter. 

It can be shown (Reference 3) that 


wr 
a(t) = [Pp(t) = Pg(t)V2y C01 + bw)" = 1] q) : 


and that for y = 1.4 and M < 2.24 this can be replaced by 


P_(t) - Po(t) 


a(t) = ——g¢Hy- (2) 
= Agee ge at i 6: 
where G(M) = 1+ aM +aqM + 7g5pM +... | (3) 
Also Pr/P. = (1 eh grt (4) 
For M< 1, PS = Py : (5) 


However, for M> 1, Po ? Pe so that to determine Pr (from the gage result 
Po) we must first find the Mach No., M. 
To do this, following Reference 3, we have 


p (G1 wey fT 
’s "| Beat i.) 


Thus, for M <1, we can determine M from Equation (4) since Py and P. are 
known (as functions of time). Then G(M) is determined from Equation (3) and 
q(t) from Equation (2). Breaking up the waveforms into time intervals (the 
same set of intervals for Fr and both Pc gages) and using this method then 
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yields a set of values of q (t,). Then Q = fq dt 3 tq(t, At, (7) 
i 


For M > 1, we can determine M from Equation (6), then Py from Equat* (4) 
since P_ and Pe are known (as functions of time). The same procedure as above 
then determines Q. (y is known as a function of the static overpressure, see 
Reference 11.) 

In practice, the overpressure-time rather than the absolute pressure- 
time histories are the data provided. Therefore, we wished to transform the 
above relations to a form suitable for using the measured pressures directly. 
To do this we rewrite Equations (4) and (6): 


PT 
P * f, (M, ¥)s M<] (4a) 
Pp 
7" fo (M, Y)> M> 1] (6a) 


and note that Equation (2) can be written: 


APy(t) = APe(t) 
q(t) = —+0) (2a) 


Equation (4a) may now be written: 


AP, + Py 
re i ee 
or, AP, ans (M, y) APo + Po [f, (Ms y) - 1], M<1 (8) 


and similarly, APy = fo (May ) APs + Py [fo (M, y) - 1], M>1. 


rar : 


Vy 


Therefore, for fixed values of M, y and P, the desired relations, Equations 
(8), are straight lines. We plotted these relations using AP» as ordinate, 
AP. as aoscissa, calculating f, and fo from the right hand sides of Equations 
(4) and (6) for assumed values of M and y. We made separate linear plots for 


ED FLA! 


each decade of static overpressure, APc. This is convenient since y can be 
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{ 
| 
obtained as a function of static overpressure (Reference 11). For 1 <P. < | 
100 psi we took y = 1.4. For 100 < Pe < 1000 psi we took y 2 1.35. 
It turns out that there were few instances of concern to us for which 
Pe > 100 psi. Further, it turns out, as will become apparent when we con- } 
sider inaccuracies involved in the data reduction process, that changes in 
both y and P, induce changes in the calculated results which are negligibly 
small compared with unavoidable inaccuracies in the data reduction process. 
We assumed P, = 14 psi. Also for M2 1.3, the difference between Po and Py 
is completely negligible compared with inaccuracies involved in data reduction. 

We plotted Equations (8) in the manner indicated for values of M between 
0.2 and 2.0 in increments of 0.2. However, in reducing the Operation Teapot 
data we only rarely encountered a value of M > 1.2, so that the Mach correc- 
tions owing to perturbation of the flow by the measuring instrument are, in 
general, quite modest’ See Equation (3). 

Prior to initiating the data reduction process we carefully examined 
the pressure-time waveforms (Reference 12) in order to devise a systematic 
method with as good accuracy as can be achieved. We noted several sources 
ef inaccuracy: the waveforms are small, especially on,the pressure axis, 
typically 0.5" peak value; the areas under the pressure-time curves are 
small, usually less than 1 square inch, this being especially true for the 
(already reduced, but not Mach corrected) dynamic pressure-time curves where 
the areas under the curves are typically in the range of 0.2 to 0.4 square 
inches (with a few smaller and a few larger than this); in many cases the 
waveforms are highly irregular, some of the irregularities probably being 
real, i.e@., physically present at the gage, and some probably being instru- 
mental; there are some problems which appeared to be due to baseline drift 
during recording and some which may well have been due to gage clogging (with 
dust). ; 

The previous data reduction process, which provided the available dynamic 
pressure-time histories, used a point-by-point method; i.e., values of total 
and static pressures were taken at identical (or almost identical) times from 
the raw data and the dynamic pressure calculated at each such time (as the 
difference between the total and stat‘c overpressures). Sometimes the dynamic 
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® In reducing the data for the Plumbbob, Smoky-15 event, Equations (4) and 
(6) were used. Plots of these Equations are given in Figures B.1 and B.2. 
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pressure, as obtained in this manner, was negative — owing presumably to 
the difficulties enumerated in the previous paragraph. 

Any attempt to use the pressure-time waveforms in the form in which they - 
are now available to obtain a dynamic pressure-time curve via a point-by- 
point data reduction methoc seemed to us to be doomed to failure. We con- 
cluded that better accurac_ (though not nearly as good as could be desired) | 
could be achieved by use of a planimeter to measure areas under the curves 
in selected time intervals. 

We used the sets of parametric curves (straight lines) corresponding to 
Equations (8) in conjunction with visual inspection of each set of waveforms 
(3 waveforms: Py and two Pos only one of the latter having been used in the 
previous data reduction process) to determine approximate average values of 
Mach No. over trial intervals. Where average values of M were relatively 
large or rapidly varying the time intervals were taken fairly small. We 
also looked visually at the waveforms in doing this in order to properly 
average out effects of irregularities in planimetering the areas under the 
curve for each time interval. When no (or only small) irregularities were i 
present we took the relative contribution of each time interval to the total 
area under the pressure-time curve into account in choosing the intervals. 
Finally, since the waveforms had to be considered in sets, if one waveform 
had much more irregularity than others in the set of three, it largely 
determined the choice of time intervals, the same set of time intervals 
being used, however, for all three waveforms in each set. The number of 
time intervals used for each set varied from ] to 6 in our data reduction 
procedure for Operation Teapot. We believe that the choices of intervals 
in each case are about the best we can make: using a greater number of 
intervals decreases the areas associated with them and results in greater 
inaccuracy in measuring the area with a planimeter. (In general, we 
planimetered each area twice, and in the case of the smaller areas three 
times — without resetting the planimeter; this improves the overall] accur- 
acy somewhat. Also, in many cases, including all those in which results 
could be questioned for any of several reasons, areas in question were 


replanimetered. ) 


B.2 PROCESSING OF REDUCED DATA 
Having obtained the areas under the various pressure~time curves for 
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each selected interval and for the total curve (value obtained separately and 
checked against the sum over the intervals) we proceeded as follows to obtain 
(4) values of dynamic pressure impulse. First we calculated average values 
of AP, and AP for each gage and for each time interval simply by dividing the 
planimeter value AP ,Aty» for example, by Aty,s the selected time interval. We 
then read the Mach No., My» for the time interval from our set of parametric. 
straight line plots of AP, versus AP. corresponding to Equations (8). Inter- 
polation between the lines of constant M, spaced at intervals of 0.2, caused 
negligible inaccuracy here compared with other inaccuracies in the data reduc- 
tion process (alluded to previously and to be discussed later). Then Equation 
(3) — or a plot thereof — determines G(M,), the Mach correction factor for 
the ith time interval. Considering now the AP. and one of the Py gages, we 
did this for each value of i and then formed (OP, - APo), dt, for each i and 
also 7 (AP, - AP<); At,. Then we associated @ weighting factor with each tine 
interval: 


apes (AP, - aPc), At; sh 
i Fab, i), at; 

This factor is the ratio of the dynamic pressure impulse associated with the 
ith time interval to the entire dynamic pressure impulse at the gage location. 
(When the ground static overpressure gage is used for AP, the two gages are 


a few feet apart; however, the ground gage does not perturb the flow field.) 


Note that tW, = 1. (10) 
In using (9) we have assumed 

AP, = APT s M2 1.3 (11) 
since, as we have previously noted, the difference between AP» and AP is 
negligible compared with unavoidable inaccuracies in data reduction. In fact 
in the Operation Teapot data reduction, use of (11) is justified in all cases. 
(Remember also that by use of time intervals rather than a point-by-point data 
reduction process we smooth out irregularities in Mj — real or instrumental 
and possibly made worse by data reduction procedure inaccuracy — and substi- 
tute an average value of M; for the it time interval; in a point-by-point 
method some higher, and lower, values of M would be found in each of our time 
intervals in which there are significant waveform irregularities. ) 
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With the aid of Equation (9) we are now in position to determine an 
average value of the Mach correction factor over the entire pressure-time 
history. Denoting this by G(M), we have 


G(M) = EW, G(M,) . (12) 


| 
| 

For reasons that will be clear when we discuss accuracy of the metiod in | 
more detail? this is the only use we made of the planimetered areas of the | 
selected intervals for each set of waveforms; in the remainder of the data 
reduction process we used only total area under the pressure-time curves. 

Note that G(M) as given by Equation (12) is greater than the value which 
would be obtained if only one “interval” — the entire waveform — were used. 
That is, if we were to base the average Mach correction factor on the areas 
under the entire P, and Ps waveforms, we would get too low a value. 

This is because the highest values of M and G(M) occur in intervals of 
time which have the greatest differences between AP and APs j.e., which 
make the greatest contribution to dynamic pressure impulse. Since the Mach 
Nos. involved are fairly low for Operation Teapot events, values of G(M) are 
also fairly low so that the difference between G(M) as given by (12) and the. 
value of G(M) which would be obtained based on the entire pressure pulse — 
one interval — never exceeds a few percent. We now can. usé the value of 
G(M) obtained from Equation (12) in two ways, separately for the two static 
overpressure gage results, each used in conjunction with the single total 
head overpressure gage result, to obtain 4 estimates of the value of dynamic 


pressure impulse. We denote these estimates by Qy° izle... 4. 
Q = f(dP, - APs.) At/G(My ) (13) 


where the integral is the difference of the planimetered values of the areas 
under the entire curves of total head overpressure and static overpressure 
(APp 2 APy for Operation Teapot events). G(M) is found from Equation (12). 
a is the value of dynamic pressure impulse associated with the static pres- 
sure gage placed at the same location as the totai head overpressure gage 
(and G(M]) is also so associated]. 

Qo is given by a similar equation, the only difference being that the 
planimetered value of AP. At and the value of G(M) are associated with the 
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ground level static overpressure gage; we use the subscript 2 to connote this. 


Then 

Qp = S(AP, - APso) at/G(M,) (14) - 
Next 

Q, = sa(t) at/G(M, ) (15) 
and Q, = fq(t) At/G(M,) : (16, 


Here fq(t) At is the planimetered value of the area under the dynamic pressure- 
time curve provided by the previous (point-by-point) data reduction. 
G(Mj) and &(Mo) rarely differ by more than a few percent in this report, 
so that Q3 and Q4 likewise rarely differ by more than (the same) few percent. 
Since all 4 estimates are dependent on a single gage result, and 2 of the 4 
estimates are based on one static overpressure gage, 2 on another, the 4 esti- 
mates are certain'’y not indapendent: however, as will soon be clear, inaccura- 
cies in obtaining dynamic pressure impulse stem not only from the raw data but 
also from the data reduction process itself. 
Thus Qy and Qo are both dependent upon the single total head overpressure 
gage result but are otherwise independent. Q3 and Gq are almost completely 
dependent upon one another, differing only because of differences between 
(My) and G(Mo), but from a data reduction standpoint are almost independent 
of Q; and Qg. Put in another way, Q; and Q3 are only dependent upon one 
another through G(M,), which has the same effect on both, and Qp and Qy are 
only dependent upon one another through G(Mo) in carrying out the data reduc- 
tion process. 
In many cases, after the estimates of dynamic pressure impulse had been 
obtained, one or more of the four estimates were either missing (gage failure, 
or missing APc gage at the ground range in question) or obviously bad. When 
a reading was obviously bad a reason for this was almost always apparent. 
Most usual reasons were: sAP_At and LAP cAt both considerably larger than the 
differences between them, which greatly magnifies errors in the data reduc- 
tion process; baseline drift or instrumental noise problems associated with 


one or more of the gages; a‘ea under curve too smal] and/or irregular ror 
accurate planimeter values to be obtained especially in the case of Sq(t)dt. 
Rather than attempt any, necessarily partly subjective, weighting of the 
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results obtained we, therefore, simply averaged the acceptable Q values (from 
1 to 4 of them) in each case. We also made a comparison of results obtained 
with values of dynamic pressure impulse inferred from displacement (based on 
data from ideal shots, Reference 1). While these data have considerable 
scatter, values of dynamic pressure impulse inferred from side-on and face-on 
dynamic pressure impulse versus displacement curves are almost always within 

a factor of 1.5 of one another. We did not use any of these results to adjust 
any values of Q; obtained from the pressure gage data. (See, however, Section 
2.2.4 for a discussion of displacement-inferred values of dynamic pressure 
impulse; here we have recourse to a larger data base than in the previous work, 
Reference 1.) 

Averaging the acceptable Q values obtained gives a somewhat higher weight 
to Q3 and Qq than to either of Q] and Qo; Q3 and Qq correspond to the previous- 
ly reduced data (point-by-point) with our Mach correction applied. As a 
comparative check of the results we plotted Qo versus Q], Q] versus Q3, and 
Qo versus Qy. Ideally the results should then lie on 45° lines. We found 
that: 

1) For Q; versus Q3 the points are scattered uniformly about a 45°. 
line; | a oo . 

2) For Qo versus Q] and Q» versus Q4 plots (with Q2 as ordinate) 
the points are scattered about a slightly greater than 45° line; 

3) The scatter of the points is about the same on all 3 plots. 
This shows that the Q; values are of equal reliability. The fact that Qo is, 
on the average, a little higher than the other Q; values may be due to a 
physically real difference resulting from the APp and APs gage locations, or 
it may be associated with the fact that the ground level APs gage does not 
’ perturb the flow whereas the flow is perturbed at the other &Ps gage location. 
Thus, we cannot say that Qo is less reliable or less likely to be correct 
than the other values of Q;- 


B.3 ERROR ANALYSIS 

Finally we must discuss the question of accuracy in more detail. We 
start with the planimeter result. When we measure the area under the entire 
pressure-time curve, a value expressed in square inches is obtained. This 
must be converted to pressure X time (psi-msec in this case although we con- 
vert the final results to kPa-sec). So in addition to the planimetering 
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error, there is an error in the scale length conversions on the raw data plots, 
psi per inch on the ordinate and msec per inch on the abscissa. We estimate 
typical measurements errors to be about 1% in the time scale, 2% in the 
(smaller) pressure scale, and ~3% in planimetering the small irregular areas, 
in spite of going around the areas more than once which diminishes part of 
the planimetering error; (it doesn't necessarily diminish the part of the 
latter error associated with inaccurate following of the contour, but it 
diminishes any planimeter resetting and vernier reading errors). In. measuring 
areas under the curve for intervals within a pressure-time waveform there is 
a further error of placing the vertical lines which subdivide the waveform 
into intervals. That is, there is an error in At, which in turn causes error 
in AP At, and APAt;, as determined with the planimeter; since At; for the 
APs gage waveform wil] not be precisely equal to the measured At; for the AP, 
waveform another small error is introduced, perhaps ~1%. We note also that 
the area within an interval is smaller than the total area under the waveform, 
usually by a factor of 2 or 3 for the most important intervals within a wave- 
form, so that the percent error in planimetering the interval areas is larger 
than for the entire area. (This is the reason for not using a larger number of 
intervals than we have used. ) , .o 

Thus before any calculations are undertaken, several small errors are: 
introduced (in addition to instrumental errors in obtaining the raw data 
with which we begin.) The root mean square value of the errors discussed is 
typically ~4% for the entire waveform and perhaps about 5% for an important 
interval within a waveform — for each of AP ,At and APcAt. The absolute error 
bounds, however, can be anywhere from zero to 6 or 8%, for the waveform and 
an interval, respectively, depending on whether the several separate errors ; 
tend to cancel or reinforce one another. 

The principal reason for the much greater inaccuracy in the final 
results achieved is due to the magnification of this error in the data reduc- 
tion process. Dynamic pressure impulse is dependent upon the difference 
between AP) At and AP.At (for an interval or the whole waveform). With the | 
estimate of typical error inputs given in the previous paragraph we can make 
a tabulation of expected and (approximate) upper bound errors as a function 
of the ratio of APsAt to APpAt. We obtain for an entire waveform: 
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% ERROR 
APsAt/AP at RMS Approximate Upper Bound 
0.25 5 10 
0.50 9 18 
0.75 20 40 


For an important interval within a waveform, the entries in the last two 
columns should be increased -y 25 to 40% (of the listed % error). For less 
important intervals, i.e., those that contribute small amounts to the total 
dynamic pressure impulse, still larger percent errors would be obtained. In 
each case the lower bound of the percent error is essentially zero, owing to 
possible cancellation of input errors. Random addition or cancellation of 
errors was undoubtedly involved in the elimination of obviously unacceptable 
results in some instances. 

We have also combined the above analysis with inaccuracy caused by our 
use of Py = Py (for M< 1.3), fixing Py = 14 psi, and interpolating between 
the parametric curves (straight lines) given by Equations (8) to determine M. 
The inaccuracies involved here are completely negligible compared with the 
unavoidable inaccuracies just discussed. Likewise errors in G(M,) and G(M) 
owing to the assumptions of convenience just mentioned are negligible compared 
with the unavoidable inaccuracies. - 

The errors we have discussed so far relate to determination of Q] and 
Qo. The errors we introduce in determining Q3 and Qq stem from two sources: 
(1) planimetering of the dynamic pressure-time waveforms (already provided 
by a previous, point~by-point data reduction process, Reference (2); (2) deter- 
mination of G(M). 

The planimetering error involved here is larger than the error in plani- 
metering AP¢ and AP, since the waveforms provided are smaller. However, there 
is no magnification of this error — this magnification has already taken 
place in forming APp - APs in the previous data reduction from which we have 
q(t) waveforms. We cannot assess the detailed sources of error involved in 
obtaining q(t). We roughly estimate, however, that errors in Q3 and Qq range 
from 7 to 10% when the area under q(t) curve is >0.5 square inch to perhaps 
30% when the area under the curve is ~0.1 square inch. From our plots, pre- 
viously mentioned, of Q2 versus Q], Qo versus Qq, and Q] versus Q3, which 
show about equal point scatter in all cases, we conclude that the inaccuracy 
in our method based on planimetering areas and the inaccuracy in the point- 
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by-point method are roughly equal. 

The fact that the errors are unavoidably as large as they are, and are 
still larger for individual intervals within a waveform, is the reason we 
used the areas under the total waveforms, Mach corrected with G(M), [with 
G(M) based on interval planimetering as previously described] to determine 
Q] and Qo — rather than using the planimetered areas AP,Aty and APgAty along 
with G(M,) to obtain the contribution of each interval to Q] and Qo, then 
summing over the intervals to obtain Q] and Qo. To show that the method we 
used was the more accurate of the two possible methods we simply explain 
why it is that we can determine G(M) accurately from planimetering values 
within intervals, even though the latter values are themselves rather inaccu- 
rate. The reasons for this are as follows: 

1) intervals with sufficient planimetered area to have high weight- 
ing factors, W;, associated with them can be planimetered much more accurately 
than small intervals; 

2) G(M) is only “large” for high values of M and these are always 
associated with relatively large planimetered values; 

3) the APpAt; and AP dt; values: will not be nearly equal for any 
interval with a large vatue of Wz [except where M is small for the entire wave- 
form, in which case G(M) is not much greater than unity, j.e., correction for 
Mach effect is small]; 

4) our weighting factor is normalized so that IW; = 1 so that if 
one or more values of W; are too high, others are necessarily too low thus 
providing some cancellation of errors; also in the important parts of a wave- 
form G(M) varies less than over the total waveform and an error in one mea~ 
sured interval is most likely to be compensated in an adjacent interval, thus 
augmenting the error cancellation effect in determining G(M). 

Had we decided to use the individual interval planimetered results, thus 
determining q(t; )dt; for each interval, some of the aforementioned cancella- 
tion would also have taken place. However, since use of the individual 
interval planimetered areas provides an accurate determination of G(M), the 
most accurate we can achieve, as previously explained — the method we have 
used yields as accurate a result as we can obtain without need for further 
use of the planimetered results for the individual intervals. 
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APPENDIX C 
(This Section is Unclassified) 


THE LEAST SQUARES PROCEDURES 


C.1 INTRODUCTION 
All of the least squares fits in this report were obtained with the 
TI-59 calculator using library programs supplied by Texas Instruments. 
All of the fits were either of the form 


y =m +b (1) 


or Y = ap + a4X) + ayXy (2) 
where y was the natural logarithm of the desired quantity; x was either the 
independent (controlled) variable or its logarithm and the least squares 
variables Xy» Xo involved only the special case xX, =x and X» = x2, 

The least squares equations are obtained by minimizing the sum of the 
squares of the deviations of the data from the curve. If the experimental 
points have a variable scatter for a given small range of x, weights should 
be applied in the procedures so that the least squares equations correspond- 
ing to Equation (1), for example, are obtained by setting the derivatives of 


S = IWiLyy - (mx, + b) 7? (3) 


with respect to m and b equal to zero, thus obtaining two linear algebraic 
equations which can be solved for m and b. Here S is the weighted sum of the 
squares of the residuals about the fitted curve and 

Wy = o*/o,” 
where o is a constant (to be determined from the deviations of the data from 
the fitted curve) and ay” is a measure’ of the expected deviation from the 
true value for an observation y; (of unit weight). 

In the Texas Instruments programs W; = 1. In individual cases to be 
discussed we will point out that the programs used are nonetheless quite 
adequate for our needs. This is mainly due to the fact that in all cases, 
we minimized the sum of the squares of the logarithmic deviations of the 
desired function. This is equivalent to minimizing the sum of the squares 
of the percentage deviations from the curve since 


A(én y) = ra (4) 
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That is, for small deviations from the curve, the fractional error dy ts 
equal to A(2n y). For large deviations this is not exact and in our tabula- 
tion of results we show both 


Jj - Yi = Yo? 


nee y; - fn y,) and ny 
Ye 


where y, = the value of y(x;) corresponding to the curve fit? 
The statistical analog of Equation (4) (Reference 13) is 


Var (2£n y) = 5 Var y (Var = variance) . (4a) 
y 


Now in all of the cases in which we fitted the logarithm of the desired 
function y, the values of y covered a large range. If we consider, for 
example, a displacement measurement of 100 metres to be in error by 2%, the 
square of its deviation from the true value is 4; if we consider a displace- 
ment measurement of 1 metre to be in error by 20%, the square of its devia- 
tion from the true value is 0.04. So even though the percent errors in our 
data are not necessarily uniform (as best these errors are known), assuming 
them to be uniform is much closer to reality than any other assumption we can_ 
make. From the example just given we readily see that the coefficients in the 
least squares fits would be determined almost entirely by the data with large 
values of y if we were to minimize absolute rather than percent deviations 
from the fitted curves. 

C.2 DYNAMIC PRESSURE IMPULSE VERSUS DISPLACEMENT (a) 

DISPLACEMENT VERSUS DYNAMIC PRESSURE IMPULSE (b) 

In application of our results we need fits of the data for cases (a) 
and (b), i.e., with each quantity used in the role of independent (controlled) 
and of dependent variable. 

In each subcase of this case the data are fitted very will by straight 
lines on log-log plots. Thus for fitting on the TI-59 calculator the fit is 


of the form y =m +b (1) 
with yzn Ig, x= fn Dd for Case (a) 
and yD, x= ln Ig for Case (b) 


* In most cases we omit the subscript i denoting the ith data point; the 
summations are taken over the data points in all cases so that omission of 
the subscript i will not lead to any confusion. 
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and b=e=2nB. (5) 
Equation (1) can then also be written [for case (a)], Iq = Bp" (la) 


with a similar relation for Case (1b). However, after we have fitted case (a) 
we can calculate (b) directly (without fitting it) provided we recall certain 
data summations from the machine memory. For Equation (1) we obtain from the 


OE a” a I CAN Rat 


TI-59 fit: 
slope = m= Exy - UxEy/N_ (6) 
2 ~ (2x)2/N 
= b = ry =. mzx (7) 
Yintercept N 
2 1/2 
and the correlation coefficient, r =m (+5) : (8) 
o 
y 


r2 is called the coefficient of determination. The various summations can be 


recalled from the machine memory. a6 and af are the variances of the x-array \ 


and y-array data and are given by * 


oy “i Z(xy - x)? (9) 
oy? =H ly; > H? (10) 
with i= x EX (11) 
and ye i ry (12) 


N is the number of data points. 
Equations (6) to (8) are not symmetrical in x and y and the algebraic in- 
verse of the least squares fit Iq = BD" (1a) 


is C= eas but this is not a least 
squares fit to the data with the role of Iq and D reversed. 
However, manipulation of Equations (6) to (i0) leads to the equation 


* Bar over a quantity indicates an average value. 
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-1/2 
r= Cay = ye Exty] (Cee? - j (2x)? aay? - f (11) (8a) 


which is symmetrical in x and y, i.e., is unchanged when x and y are inter- 
changed. 

Thus, if we have obtained the least squares fit Iq ® BD" we know that 
the correlation ratio r has the same value with the roles of Iq and D 
reversed. It can be shown (by taking the origin at the data centroid), (Refer- 
ences 13 and 14), that the least squares fits 


= ann 
I, = BD (1a) 
meee 2 (1b) 
D Al, 
have values of 2 and m which satisfy Lm re : (13) 


Thus, having determined m (and r) in the least squares fit to case la 
we obtain 2 in the least squares fit to case (1b) from Equation (13). The 
constant A in Equation 1b is then determined by the fact that both least 
squares fits (1a) and (1b) pass through the centroid of the data distribution. 
(Compare Reference 15.) Thus, from Equation (1b) we have 

mdz hHA+L&n Ig » * (1b) 

and in particular Ir D= MKMA+£2n I, (14) 
where from Equations (11) and (12) 


en 


So 
a 


Een D) (15) 


Zi Tq =} 2(en Iq) 


and these values are available from the machine memory after running case 
(la). Thus from Equation (14) we determine 2n A, hence A, and the least 
squares fit to case (1b) is fully determines. 

Having obtained the least squares fits, the coefficient of determination 
provides a measure of the goodness of fit. By manipulation of Equation (3), 
with W; = 1, and Equations (6) through (12) it can be shown that 

T-r sy, (16) 
Noy, 


In this report we are concerned only with positive correlation, so that the 
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nearer r is to +1, the smaller is S, the sum of the squares of the deviations 
from the fitted curve, and hence the better the fit to the data. The expres- 
sion for 1 - r@ can also be put in the form 


l-rs - (17) 
thus expressing 1-r? as a ratio of the sum of the squares of the deviations 
from the curve to the sum of the squares of the deviations from the mean (of 
the data). The denominator may be regarded as a normalizing factor. It is 
independent of the functional form used for the fitting function. Without 
such a denominator 1-r? would tend to increase with the number of data points 
— even if the data were excellent. 

Now letting Son I and Son D be the sum of the squares of the deviations 
from the fitted curve when 2n Ig and én D, in turn, play the role of dependent 
variable, we have from Equation (16) 


S S 
(1 - r?) = “en Iq = £n_D (18) 


2 
. Non tg end 


Here a ig n HE (in ig)? - i (ren Tad 
(19) 
oy. ale WE (en p)? - r (en D)*3 


in accord with Equations (9) through (12). Since we have seen that r is 
unchanged when the roles of D and Ig are reversed, Equation (18) shows that 
the sums of the squares of the deviation from the curves are not the same for 
the two ccrresponding least squares fits (to the same data) but that 


Sent oon I 
- + = 2g (20) 
fn D Cen D 


In all cases treated in this report we found Sen iF to be considerably 


less than Sey D’ 
In an application to be made of the results of this case we wish to 
know the standard deviation of °’ ie value of the slope of the curve as well as 
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the standard deviation of our observation. We use the following notation 
with the subscript c referring to the curve fit in each instance: 


on test = variance of an observation of én Iq relative to the 
curve; (note that we have assumed the variance of 
the percent error to be independent of the abscissa 
value) 


a, of = variances of slopes of curve fits. 


Then it can be shown that (Reference 16) 


2 Zz 
Sen IgsC * Sen Iq/ N-2* 


j (21) 
Sen D.C Sen »/ N-2 
5 -2 2 9 
ot = fit D . ir fn D 
, . £L No; : N- of 
en Ig en Iq 
(22) 
2 
S -2 20 
ot = en A a ler" “kn I 
m N-2 
Nop, D of 
fn 0 
Using Equation (8) twice, first with D, Iq being independent and dependent 
variable and then with their roles reversed, we also see that 
o 
£2 
—2= : (23) | 
a, =m 


In the tabular results for the Displacement - Dynamic Pressure Impulse 
data curve fits, we list for each subcase the quantities N, Sen la’? [en iade: 
Spy D? en D,C? 2, m, Tos Oy in addition to r and the fitting functions. 


To gain further insight into the reliability of the data fits we have 


* Here there are N-2 degrees of freedom; tne two degrees of freedom lost 
correspond to the number of regression coefficients (two). 
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also tabulated the following additional quantities defined as follows: 


D-Dc ,2 
Si. = i(——= 
I,-I aa 
q-*4c 
S$; =I 
Ig ( Iq ) 
Sh 1/2 
9p,C (q=2) 
25 
flay _ 
Ot = 
Iq ~ Sep) 
and finally the root mean square percent error which is 
Ey = (Sin)'/2 x r00% 
(26) 


E, = (st sn)'/? x 1002 
Tq Ty 


The o' quantities, Equation (25), are standard deviations of the fractional 
error (percent error apart from a factor of 100). Thus when the number of 
data points N is fairly large (large enough so that the fit is fairly reli- 
able) the values of Op» Of are aimost the same as the values of Ep and EIg? 
respectively, (aside from the factor of 100 expressing the latter quantities 
in percent). 

The reason we have chosen to list these additional quantities, in this 
and in other data fits to be discussed shortly, is that in many cases the 
data scatter is rather large (irrespective of the functional form selected * 
for fitting the data). This means that a deviation Afn y in £n y for a given 
data point may differ considerably from Ay/y. To illustrate, suppose the 
curve fit and data point values of a displacement are 10 metres and 6 metres. 

The contribution of this point to the sum of the squares of the deviations | 
of £n D is then (fn 10-£n 6)2 = 0.261 (irrespective of which value is the 


rea curve fit values and which is the data point value). 
seine The contribution to S) is 9 

ae (1056)" = 0.160 
atte mI) : 
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if the curve fit value is 10 and the data point value is 6; if, however, the 
curve fit value is 6 and the data point value is 10, the contribution to Sp) 
= (10-6)" = 0.444. 

Since the least squares curve tends to pass through the data region with 
a fairly uniform distribution of points on either side, on the average the 
curve fit will be below the data point about as often as above it. Averaging 
the above results yields 0.160 0.7 = 0.302 which is to be compared with 
the value 0.261 for the logarithmic deviation. Had we chosen the values 10 
and 9.8 instead of 10 and 6, the two compared numbers would be virtually 
identical. 

It is inherent in the nature of the data we are dealing with in this 
report, as is clear from the discussion of the data, that there is a fairly 
large data scatter in many cases. Also in many cases the number of data points 
jis rather fewer than we desire. In cases where the scatter of the data is 
smal] we find that Sa is fairly close in value to Sp, p and similarly S;, is 
fairly close to Sp, Iq° Usually Sen p< Sp» as in the case illustrated, and 
Son os Sy} ; the situation is sometimes reversed, however, usually because 
one or two data points with large deviations from the curve fit are below the 
curve. Conversely, for a case in which one or two data points are far above 
the curve it can happen that Sy >> Son y since there is no limit to the con- 
tribution such a point can make to S', y here indicating any quantity whose 
gn has been fitted. (For a data point helow the curve fit the maximum con- 
tribution to Sy is 2 2 


Finally, quantities such as OF Cc and Er are directly associated with 
the plotted data and are easily visualized. Since they pertain to fractional 
or percent errors they apply equally well at all points along the curve fit 
(although the curve fit has greater predictive accuracy in the vicinity of 
the data centroid than toward the ends of the data range). On the other 
hand Cnn I refers to the natural logarithm of the plotted data (ordinate) 
and the latter may be quite small or large (or even negative), so that while 
py ia is a standard measure of the data scatter with respect to the fitted 
curve, it is not as simple to interpret as are the additional quantities 
which we have tabulated. 
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In sum, the tabulated measures of goodness of fit provide some insight 
into the assessment of the reliability and utility of the data and the 
curve fits and help with the interpretation of results. 

The quantities described here are also tabulated for various curve fits 
of data other than displacement - dynamic pressure impulse data. We discuss 
the least squares treatment of the remaining data in the following sections. 


C.3 ADDITIONAL LEAST SQUARES FITS 

Displacement versus Scaled Ground Range (a) 

Dynamic Pressure Impulse versus Scaled Ground Range (b) 

Scaled Dynamic Pressure Impulse versus Scaled Ground Range (c) 
In each case the fitted curve is described by 


{ny = a +a, + aox” (27) 


where x = scaled ground range. For cases (a) and (b) several scaling factors 
were used and an approximate optimum obtained. 

In running these cases on the TI-59 calculator we feed in the natural 
logarithm of the y coordinate of each data point. Again, the machine gene- 
rates a fit and all of the relevant summations over the data are available 
from the machine memory. For this type fit (trivariate), however, the 
machine provides a quantity R© rather than r as a measure of the goodness of 
fit. Ris the multiple linear correlation coefficient between y and the other 
least squares variables xX, =X and Xo = x*; Rn? is called the coefficient of 
determination (see Reference 16) and is given by 


2 2 
ro, try -2%,5 fr 7 
Re = yl y2 a yl _y2 (28) 
ea 
The subscripts 1 and 2 refer to x] and xa; ry ry2 and Ty are given by 
equations analagous to Equation (8a). a 
2 2\ 
ry 7 Lexy ~ RBA tax? - GN" aay? - 3) 
2. rxés A fixe hans 2 (ty) hg a 
ry = [ex’y - aca (ex Pepe Say i (29) 
be 
2 2 2,2 \-1/2 ia 
r = frx? . UXIX ] [rx? _ (=x) Itzx" BS (=x ) 1] i 
12 N N N 
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(The superficial lack of symmetry between "v1 and rye results from the fact 
that the least squares variables x] and xg as used here are x and x2, respec- 
tively.) The quantities rr Tyo and Pig are the coefficients of correlation 
(also called simple correlation or zero-order coefficients) between y and x7, 
y and x9, and x] and xg, respectively. The important point here is that 
2 >» 

(1 - R°) = (17a) 
so that R2 provides a measure of the goodness of fit precisely similar to that 
provided by r? in the previous discussion: Re is a measure of the closeness 
of fit of the regression plane (in £n y, x], x2 space) to the data points. 


There is one further change: replacement of N-2 by N-3 in Equations (21) 
and (25) so that 


2 = - 
St 1qi6- Sen qh 3 


2 . (21a) 
Cen Dd, C Sen pfs 


and 
| ope = (Sp/t-3)' 
1/2 a 
+= (St INS 
ae Sige 


We list the same quantities in the tabular results rel~.ting to the 
fitted curves in this subsection as in the previous subsection. 

Some additional fits were required in the cases treated in this subsec- 
tion owing to the form of Equation (27). This functional form entails either 
a maximum or a minimum (which may fal] within the data range, near one end 
of the data range, or at a point remote from the range of the data). In the 
case of Dynamic Pressure Impulse versus Scaled Ground Range and Scaled Dynamic 
Pressure Impulse versus Scaled Ground Range, the experimental data are such 
that the curves are concave upward and have a minimum usually approximately 
at the large x (= scaled ground range) end of the data range. In these 
instances, in order to have a fit of the proper shape for large x, and capable 
of extrapolation with no more than the usual uncertainty in extrapolating a 
monotonic well-behaved function somewhat beyond the range validated by data, 
we generated additional fits (of the same functional form and with the same 
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scaling in each instance) over the data contained in the large x end of the 
data range. Thus, in these cases we have two fits of the form of Equation 
(27), one over the entire data base and one over the portion of the data base 
relevant to large scaled range. 

The need for the additional fits stems from the fact that there are many 
targets susceptible to damage at lower levels of dynamic pressure impulse than 
are the vehicles of interest in this report. 

There are possibly other methods of accomplishing this purpose, includ- 
ing use of functions other than provided by Equation (27) for fitting the 
data. We note, however, that the data itself is concave upward, so that any 
function conforming to this shape and consisting of more than two terms all 
either increasing or decreasing monotonically would seemingly have at least 
one positive and one negative term and therefore a possibility of a minimum 
somewhere (though not necessarily near or within the range of the data). In 
an attempt to avoid this we used a one term expression for the fit in one 
case: 

y = a exp[-b(x-xq)"] 

We chose a and xg so that the curve was forced to pass through a selected 
point: The exponent m allows some flexibility in adjusting the curvature as 
a function of x as compared with a simple exponential and y decreases mono- 
tonically as x increases for x > xg. (x = Xg was chosen as the lowest value 
of x in the data range.) Then with a change of variable Y = en(2) and X = 
X=Xq the fitting equation takes the form fn Y= Mm b+mnX. We fitted 
this equation to the data for one case. It did not fit well at all. (At 
least part of the reason for this, however, is that the fitting equation is 
ill-conditioned to fitting the data without use of a weighting function — 
large percentage errors in y for smal] y have little effect on the data; as 


previously pointed out the Texas Instruments library programs for least squares 


fitting do not allow use of weighting functions. ) 

Thus, while it is possible that some other (probably not more than two 
term — however complicated) function or some transformation of variables may 
allow good fit to the data without a minimum (or maximum) in or near the data 


range we believe it simpler and quite justifiable to use Equation (27) for all 
of the fits in this subsection in the manner we have described. Equation (27) 


is a simple and convenient function to use. 
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APPENDIX D 
(This Appendix Is Unclassified) 


LIST OF ABBREVIATIONS, ACRONYMS, AND/OR SYMBOLS 


Ip - static pressure impulse 
Iq - dynamic pressure impulse 


Ig, - dynamic pressure impulse: initial estimates for ideal /near-ideal 
events; first iteration results for non-ideal events 


Tq - final estimates of dynamic pressure impulse: iterated results for 
ideal/near-ideal events; second iteration results for non-ideal 
events 


P, - peak overpressure 
D - displacement 

da - damage 

W, Y - weapon yield 

R - ground range 

X - scaled ground range 
HOB - height of burst 
Py - ambient pressure 


Ty - ambient temperature 

Sq> Sq> Si» Sp» St - scaling factors 

n ~ exponent in scaling law 

a9» aj, a9, b, 2, m, XQ, A, B - constants in least squares fits 
x, y - independent and dependent variaples in least squares fits 
f, f1, fo, G - functional forms 

a - standard deviation 

o@ ~ variance 

4 N - number of data points 

r - correlation coefficient 

R - multiple -linear correlation coefficient 

r2, R2 - coefficient of determination 

S - sum of squares of deviations of data points from fitted curve 


S' - sum of squares of fractional (relative) deviations of data 
points from fitted curve 


x’ 
4 
> 
mY 
nS 
"4 


1/2 
S 7 at 
on See rmerer regression coefficients’ Standard deviation 


of fractional error 


1/2 
E - ) x 100 = root mean square percent error 
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- mass of vehicle 

- moment of inertia of vehicle (about designated axis) 

- projected area of vehicle 

- coefficient of friction 

0 ~ velocity, initial velocity 

force 

angular velocity 

- angular acceleration 

- torque (vehicle overturning) 

~ distance travelled by vehicle (in equations of motion) 
F(V), g(V) - functions of velocity (in equations of motion) 
C - constant in curve fits 

K, Kj,» Kp - constants in equations of motion 


x Ak Ee Mae rt rH = 
< 


Ps - absolute free stream static pressure 
Pr - absolute free stream total pressure 
Pp - absolute total head Pitot pressure 
P. - ambient pressure 


> 

~v 
nn 
fl 


Ps - Po — free stream static overpressure 
AP; = PT - Po — free stream total overpressure 
APy = Pp - Po — total head overpressure 

q - dynamic pressure 


Qj - dynamic pressure impulse; subscript ji denotes ith method of 
obtaining result (4 values of i used) 


y - ratio of specific heat at constant pressure to specific heat at 
constant volume for air 


M - local free stream Mach No. of flow behind blast front 
Aty - ith time interval (in a pressure-time waveform) 
G(M) - Mach correction factor 

G(M) - average value of G(M) overpressure-time history 

W; - weighting factor for ith time interval 
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Defense Special Weapons Agency 
6801 Telegraph Road 
Alexandria, Virginia 22310-3398 


TRC 13 April 1998 


MEMORANDUM TO DEFENSE TECHNICAL INFORMATION CENTER 
ATTN: OCQ/MR BILL BUSH 


The following reports have been reviewed by the Defense 
Special Weapons Agency Security Office: 


DNA-5056F, AD- rere YDTL- 80, 0808 
DNA-5826F, AD-C040572, "€DTL-87, 0355 C FED 
DNA-5826F-SUP, AD- C041417,“DTL-871167 CLD 


The Security Office has declassified all of the listed 
reports. Further, distribution statement "A" applies to all of 
the reports. 
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ARDITH JARRETT 
Chief, Technical Resource Center. 


